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Abstract. There has recently been growing evidence for the existence of neutron 
stars possessing magnetic fields with strengths that exceed the quantum critical field 
strength of 4.4 x 10^^ G, at which the cyclotron energy equals the electron rest mass. 
Such evidence has been provided by new discoveries of radio pulsars having very high 
spin-down rates and by observations of bursting gamma-ray sources termed magnetars. 
This article will discuss the exotic physics of this high-field regime, where a new array 
of processes becomes possible and even dominant, and where familiar processes acquire 
unusual properties. We review the physical processes that are important in neutron 
star interiors and magnetospheres, including the behavior of free particles, atoms, 
molecules, plasma and condensed matter in strong magnetic fields, photon propagation 
in magnetized plasmas, free-particle radiative processes, the physics of neutron star 
interiors, and field evolution and decay mechanisms. Application of such processes 
in astrophysical source models, including rotation-powered pulsars, soft gamma-ray 
repeaters, anomalous X-ray pulsars and accreting X-ray pulsars will also be discussed. 
Throughout this review, we will highlight the observational signatures of high magnetic 
field processes, as well as the theoretical issues that remain to be understood. 



1. Introduction 



Since their theoretical conception by Baade & Zwicky (1934) neutron stars have been 
fascinating celestial objects, both for study of their exotic environments and for their 
important place in stellar evolution. Among the first signals to be detected from 
neutron stars came from radio pulsars (Hewish et al. 1968), spinning many times per 
second, distinguishing themselves from the background of interplanetary scintillation 
signals by their extremely regular pulsations. Pulsars were also soon discovered to be 
spinning down, their periods increasing also very regularly. The rotating magnetic- 
dipole model (Pacini 1967, Gold 1968, Ostriker Sz Gunn 1969), in which the pulsar 
loses rotational energy through magnetic dipole radiation, was dramatically confirmed 
with the discovery that the spin-down power predicted for the Crab pulsar was a nearly 
perfect energetic match with the radiation of its synchrotron nebula. The rotating dipole 
model also accounts for the observed rate of spin-down of the Crab and other pulsars. 
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with required surface magnetic fields in the range of 10^^ — 10^^ Gauss for the first 
detected pulsars. This range has since significantly broadened, first with the discovery 
of a class of pulsars having periods of several milliseconds (Backer et al. 1982), believed 
to have been spun- up by accretion torques of a binary companion (Alpar et al. 1982), 
and much lower surface magnetic fields in the range of 10^ — 10^° Gauss. Recent surveys 
have also discovered pulsars with very high period derivatives (e.g. Morris et al. 2002, 
McLaughlin et al. 2003) that imply surface fields up to around 10^^ Gauss. 

Another class of neutron stars was discovered at X-ray and 7-ray energies and 
may possess even stronger surface magnetic fields. Such stars are now referred to 
as magnetars (Duncan & Thompson 1992), since they most probably derive their 
power from their magnetic fields rather than from spin-down energy loss (see Woods & 
Thompson 2005). Within the magnetar class there are two types of sources that were 
originally thought to be very different objects, although they are now believed to be 
closely related. The Anomalous X-Ray Pulsars (AXPs) were discovered as pulsating 
X-ray sources in the early 1980s and were thought at first to be an unusual type of 
accreting X-ray pulsar, from which the name is derived. The AXPs have periods in a 
relatively narrow range of 5 - 11 s, are observed to be spinning down with large period 
derivatives (Vasisht & Gotthelf 1997) and have no detectable companions or accretion 
disks that would be required to support the accretion hypothesis. Interpretation of 
their period derivatives as magnetic dipole spin down imply magnetic fields in the range 
10^^ — 10^^ Gauss. But such high fields were not widely accepted initially, since their 
detected X-ray luminosities of around 10'^^ erg s~^ exceed their spin-down luminosities by 
several orders of magnitude. It was only by connection to another subclass of magnetar, 
the Soft Gamma- Ray Repeaters (SGRs), that the extremely high magnetic fields of 
AXPs were adopted. SGRs were discovered as transient sources that undergo repeated 
soft 7-ray bursts, usually in widely separated episodes. They undergo both repeated 
smaller bursting of subsecond duration as well as much more luminous superflares, 
lasting hundreds of seconds, which so far have not repeated in any single source but 
may be repeating on much longer timescales. It was not until some twenty years after 
their discovery that their quiescent X-ray pulsations were detected and also very high 
period derivatives (Kouveliotou et al. 1998), both with a range of values very similar 
to those of AXPs. Recently, bursts resembling the smaller bursts of SGRs were seen 
from several AXPs (Kaspi et al. 2003), making it likely that SGRs and AXPs are 
two variations of the same type of object (Thompson & Duncan 1993), very strongly 
magnetized, isolated neutron stars possibly powered by magnetic field decay. 

The periods and period derivatives of the various types of isolated pulsars are shown 
in the P-P diagram of Figure ^ Assuming that the spindown torque is due to magnetic 
dipole radiation, two quantities can be defined from the measured P and P for each 
pulsar: (1) Characteristic age P/{2P): From Q oc —fl^ (where Q = 2tt/P), the age of 
the pulsar is found to be T = (P/2P)[1 — (P/Pj)^], where Pj is the pulsar's initial spin 
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period. (2) Surface dipole magnetic field: 

wfiere P15 = P/(10^^^ ss^^) (and P is in units of second), and and / (^ 10^^ g cm^ 
and R (~ 10^ cm) are the neutron star moment of inertia and radius. There are 
presently around 1600 spin-powered radio pulsars known, with periods from 1.5 ms -8s 
(Manchester 2004)|. Some fraction of these pulse at other wavelengths, including about 
10 in 7-rays and about 30 in X-rays. The magnetars, eight AXPs and four SGRs§, 
occupy the upper right-hand corner of the diagram and curiously overlap somewhat 
with the region occupied by the high-field radio pulsars. However, the two types of 
objects display very different observational behavior. The high-field radio pulsars have 
very weak or non-detectable X-ray emission and do not burst (e.g. Kaspi & McLaughlin 
2005), while the magnetars have no detectable radio pulsations, with the exception of 
the recent detection of radio pulsations in the transient AXP XTE J1810-197 (Camilo 
et al. 2006). The intrinsic property that actually distinguishes magnetars from radio 
pulsars is presently not understood. 

A third class of strongly magnetized neutron stars are the accreting X-ray 
pulsars (Parmar 1994). These sources are members of binary systems with high- 
mass companions that either have strong stellar winds or overflow their Roche lobes, 
transferring material to the neutron star. Inside an Alfven surface where the pressure of 
the accretion flow equals the neutron star magnetic field pressure, the accreting material 
is funneled along the magnetic field to the poles. The heated accretion flow then radiates 
from hot spots that rotate with the star, thereby producing pulsations. However, since 
accretion torques dominate the period derivative evolution, the surface magnetic fields 
of X-ray pulsars cannot be determined using the rotating dipole model, as they are for 
the rotation-powered pulsars and magnetars. Instead, magnetic fields of these objects 
have been measured from the energies of the cyclotron resonant scattering features 
(CRSFs) that appear in their spectra (see Orlandini & dal Fiume 2001, for review), 
since the fundamental occurs at the electron cyclotron energy, Ece = heB/{mec) = 
11.58 (P/10^^ G) keV. Table 1 lists the X-ray pulsars, the energies of CRSFs that have 
been detected in their spectra and the inferred magnetic field strengths. In most cases, 
the measured magnetic fields range from 1 — 5 x 10^^ Gauss, with the highest fields being 
around 4 x lO^^Gauss. The formation of such line features will be discussed in section 

Clearly, there is a broad range of astrophysical sources in which the magnetic 
fields approach and exceed the quantum critical field strength, Bq = m^c^/ {eh) = 
4.414 X 10^^ G, at which the cyclotron energy equals the electron rest mass. In 
this regime, the magnetic field profoundly affects physical processes and introduces 
additional processes that do not take place in field-free environments. From a physics 
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Figure 1. Plot of period vs. period derivative for the presently known rotation- 
powered pulsars and magnetars. Lines of constant characteristic age, P/2P, and 
surface dipole field (see Eqn ^) are superposed. 
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Source 



Energy (keV)'' B (lO^^ Gauss)^ 



V 0332+53 
Cep X-4 



4U 0115+63 
4U 1907+09 
4U 1538-52 
Vela X-1 



A 0535+26 



Cen X-3 
X Per 

XTE J1946+274 



MX0656-072 
4U 1626-67 
GX 301-2 
Her X-1 



12 (4) 1.0 

18 (1) 1.6 

20 1.7 

25 (1) 2.2 

27 (2) 2.3 

28 2.4 
28.5 2.5 

29 2.5 
36 3.1 

36 3.1 

37 3.2 
37 3.2 
41 3.5 
47'^ (2) 4.1 



" Data is from Heindl et al. (2004). 

Numbers in parentheses are the number of cyclotron harmonics detected. 
^ Magnetic field strength for viewing angle along field direction. 

Kretschmar et al. (2005) 

point of view, strongly magnetized neutron stars provide the only environment in which 
to measure and test these effects. Prom an astrophysics point of view, it is necessary 
to investigate the physics of strong magnetic fields in order to effectively model and 
understand the nature of the sources. This article will attempt to address both points 
of view by providing a review of the basic physical processes important in neutron star 
interiors and magnetospheres, as well as a review of the source models in which these 
processes play a critical role. 

We begin with a description of a free electron in a magnetic field in both classical and 
quantum regimes. Subsequent sections then discuss and review current understanding 
of the behavior of matter, atoms, molecules and plasma, in strong magnetic fields, 
photon propagation in magnetized plasmas, free-particle radiative processes, the physics 
of neutron star interiors, and field evolution and decay mechanisms. Then we will 
review models for magnetized atmospheres, non-thermal radiation from rotation- 
powered pulsars, burst and quiescent radiation from SGRs and AXPs, and emission 
from accreting X-ray pulsars. Other useful books and reviews on strongly magnetized 
neutron stars include Meszaros (1992), Duncan (2000), Lai (2001) and Harding (2003). 
A complimentary work concentrating on stellar magnetism in general is the book by 
Mestel (1999). 
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2. Electrons in Strong Magnetic Fields 

The quantum mechanics of a charged particle in a magnetic field is presented in many 
texts (e.g., Landau & Lifshitz 1977; Sokolov & Ternov 1968; Meszaros 1992). Here we 
summarize the basics needed for our later discussion. 

Consider first the nonrelativistic motion of a particle (charge Cj and mass mj) in 
a uniform magnetic field B (assumed to be along the 2;-axis). In classical physics, the 
particle gyrates in a circular orbit with radius and (angular) frequency given by 

micvi_ \ei\B . 

\ei\B niiC 

where v± is the velocity perpendicular to the magnetic field. In non-relativistic quantum 
mechanics, the kinetic energy of the transverse motion is quantized in Landau levels 

Ei_ = ^m^vi = (ul + ^) huj„ = 0, 1, 2, • • • (3) 

where 11 = p — (ej/c)A = m^v is the mechanical momentum, p = —ikV is the canonical 
momentum, and A is the vector potential of the magnetic field. 

For an electron (m^ — me, Cj — >■ — e), the basic energy quantum is the cyclotron 
energy 

eB 

Ece = huJce = h = 11.577 5i2 keV, (4) 

rUeC 

where B12 = -8/(10^^ G) is the magnetic field strength in units of 10^^ Gauss. Including 
the kinetic energy associated with the z-momentum [p^) and the spin energy, = 
eh/{2mf.c)(T ■ B = huJceCzf^ (with cr^ = ±1), the total electron energy can be written as 

En = nhuce + (5) 

where the index n = + {1 + <7z)/2 = 0, 1, ■ ■ ■. Clearly, for the ground Landau level 
(n = 0), the spin degeneracy is one (a^ = — 1); for excited levels, the spin degeneracy is 
two (o"2 = ±1). 

Apart from the spin degeneracy, the Landau level associated with Eq. Q is 
degenerate by itself, refiecting the fact that the energy is independent of the location of 
the guiding center of the gyration. To count the degeneracy, it is useful to define the 
pseudomomentum (or the generalized momentum) 

K = n + (ei/c)B X r. (6) 

That K is a constant of motion (i.e., it commutes with the Hamiltonian) can be easily 
seen from the classical equation of motion for the particle, dYl/dt = {ei/c){dr/dt) x B. 
The parallel component is simply the linear momentum p^, while the perpendicular 
component is related to the position vector Rc of the guiding center by 
^ cK I X B c „ - 

R. = ^^ = -^n,xB + r,. (7) 
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(B is the unit vector along B). From Eqs. and ((Tj), we see that classical radius of 
gyration, p = |r_L — (see eq. |2]), is quantized according to 



\r± - Rcl = ^in^l {2nL + l)'/'p, (8) 



where 



is the cyclotron radius (or the magnetic length). Since the two components of do not 
commute, [Krc,Ky\ = —ih{ei/c)B, only one of the components can be diagonalized for 
stationary states. This means that the guiding center of the particle can not be specified. 
If we use Kx to classify the states, then the wavefunction has form e'^^'^'^/^(f){y) (Landau 
& Lifshitz 1977), where the function 0(?/) is centered at yc = —cK^/^eiB) [see Eq. (jZj)]. 
The Landau degeneracy in an area Ag = L"^ is thus given by 

(10) 

where we have used K^^g = —CiBLg/c. Alternatively, if we choose to diagonalize 
Kj_ = + -^^5 obtain the Landau wavefunction VF„m(r±) in cylindrical coordinates 
(Landau & Lifshitz 1977), where m is the "orbital" quantum number (denoted by s or 
—s in some references). For the ground Landau level, this is (for = — e) 

WomM ^ WUpA) = (2^)1/2^ (^) -M-^m<P), (11) 

where the normalization / d'^T± {Wml"^ = 1 is adopted. The (transverse) distance of the 
particle's guiding center from the origin of the coordinates is given by 

|R,| ^p™ = (2m + l)i/2p, m = 0,1,2,- •• (12) 

The corresponding value of K± is Kj_ = {h\ei\B /c){2m + 1). Note that Kj^ assumes 
discrete values since m is required to be an integer in order for the wavefunction to be 
single-valued. The degeneracy of the Landau level in an area Ag = irRg is then 
determined by pmg — {^mgY^'^p = Rg, which again yields rUg = Ag\ei\B / (he) as in 
Eq. (fTUI) . Note that despite the similarity between Eqs. (fT^ and (jHI), their physical 
meanings are quite different: the circle p = Pm does not correspond to any gyro-motion 
of the particle, and the energy is independent of m. Also note that Kj_ is related to 
the z-projection of angular momentum J^, as is evident from the e"*'"'^ factor in the 
cylindrical wavefunction [Eq. (fTT|) ]. In general, we have 

J, = xPy - yP, = 7r^(Ki - nl) = (m - nz.)^, (13) 

where we have used Yl\ = {h\ei\/c)B{2nL + 1). 

For extremely strong magnetic fields such that huce ^ mec"^, or 

2 S 

B^Bq = = AAIA X 10^=^ G, (14) 

ea 



Physics of Strongly Magnetized Neutron Stars 



8 



the transverse motion of the electron becomes relativistic. The transverse momentum 
is again quantized according to Eq. Q, i.e., 11^ — > (2^^ + l)TieB/c. By solving the 
Dirac Equation in a homogeneous magnetic field, we find that Eq. (jS)) for the energy of 
a free electron should be replaced by (e.g., Johnson & Lippmann 1949) 



The shape of the Landau wavefunction in the relativistic theory is the same as in the 
nonrelativistic theory, as seen from the fact that p is independent of the particle mass 
(see Sokolov & Ternov 1968). Higher order corrections in to Eq. (|T5|l have the 
form {a/4:7i)meC^F{B/BQ), with a = 1/37 the fine structure constant, F{P) = —j3, 
(3 = B/Bq, for /3 < 1 and F{(3) = [ln(2/3) - {-fE + 3/2)]^ + ■ ■ ■ for /3 > 1, where 
■jE = 0.5772 is Euler's constant (Schwinger 1988); these corrections are negligible. 

Calculations based on electron wavefunctions that are solutions to the Dirac 
Equation in Cartesian and in cylindrical coordinates have both appeared in the 
literature. As is true for all quantum processes, the rates and cross sections depend on 
the choice of electron wavefunctions in a uniform magnetic field, which are dependent 
on spin state. The two most widely used relativistic wavefunctions are those of Johnson 
and Lippman (1949) and Sokolov & Ternov (ST) (1968). The Johnson and Lippman 
(JL) wavefunctions are derived in Cartesian coordinates and are eigenstates of the 
kinetic momentum operator tt = p + eA/c. The Sokolov & Ternov (ST) wavefunctions, 
given for the ground state in Eqn (|TT|l . are derived in cylindrical coordinates and are 
eigenstates of the field parallel component, fiz, of the magnetic moment operator. Given 
the different spin dependence of the ST and JL eigenstates, one must use caution 
in making the appropriate choice when treating spin-dependent processes. Herold et 
al. (1982) and Melrose and Parle (1983) have noted that the ST eigenstates have 
desirable properties that the JL do not possess, such as being eigenfunctions of the 
Hamiltonian including radiation corrections, having symmetry between positron and 
electron states, and diagonalization of the self-energy shift operator. As found by 
Graziani (1993), the ST wavefunctions also diagonalize the Landau-Dirac operator and 
are the physically correct choices for spin-dependent treatments and in incorporating 
widths in the scattering cross section. Although the spin-averaged ST and JT cyclotron 
decay rates are equal, for example, the spin-dependent decay rates are not, except in 
the special case in which the initial momentum of the electron parallel to the magnetic 
field vanishes. An additional desirable property of the ST states is that spin-dependent 
states are preserved under Lorentz transformations while Lorentz boosts mix spin states 
of the JL wavefunctions. 

3. Matter in Strong Magnetic Fields 

When studying matter in magnetic fields, the natural (atomic) unit for the field strength, 
Bq, is set by equating the electron cyclotron energy huce to the characteristic atomic 




(15) 
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energy e^/ao = 2 x 13.6 eV (where is the Bohr radius), or equivalently by p = oq. 
Thus it is convenient to define a dimensionless magnetic field strength b via 

b=^; 5o = = a^^Q = 2.3505 X 10^ G (16) 

{a = 1/137 is the fine structure constant). For b ^ 1, the cyclotron energy huce is much 
larger than the typical Coulomb energy, so that the properties of atoms, molecules and 
condensed matter are qualitatively changed by the magnetic field. In such a strong field 
regime, the usual perturbative treatment of the magnetic effects (e.g., Zeeman splitting 
of atomic energy levels) does not apply (see Garstang 1977 for a review of atomic physics 
at b <, 1). Instead, the Coulomb forces act as a perturbation to the magnetic forces, and 
the electrons in an atom settle into the ground Landau level. Because of the extreme 
confinement (p <^ oq) of the electrons in the transverse direction (perpendicular to the 
field), the Coulomb force becomes much more effective in binding the electrons along 
the magnetic field direction. The atom attains a cylindrical structure. Moreover, it is 
possible for these elongated atoms to form molecular chains by covalent bonding along 
the field direction. Interactions between the linear chains can then lead to the formation 
of three-dimensional condensates. 

Note that when studying bound states (atoms, molecules and condensed matter 
near zero pressure) in strong magnetic fields, it is adequate to use nonrelativistic 
quantum mechanics, even for B > Bq. The nonrelativistic treatment of bound states 
is valid for two reasons: (i) For electrons in the ground Landau level, eq. (fT3j) for the 
free electron energy reduces to ~ mgC^ +p^/(2me) for p^c <^ mgC^; the electron 
remains nonrelativistic in the 2;-direction (along the field axis) as long as the binding 
energy Eb is much less than nieC^; (ii) As mentioned before the shape of the 

Landau wavefunction in the relativistic theory is the same as in the nonrelativistic 
theory. Therefore, as long as ^'^/(meC^) <^ 1, the relativistic effect on bound states is 
a small correction. For bulk matter under pressure, the relativistic correction becomes 
increasingly important as density increases (see 

Note that unless specified, the expressions in this section will be in atomic units 
(a.u.), where mass and length are expressed in units of the electron mass me and the 
Bohr radius Oq = 0.529 x 10~^ cm, energy in units of 2 Ryd = e^/oo = 2 x 13.6 eV, field 
strength in units of Bo [Eq. ()16|l ]. A previous (more detailed) review on this subject is 
Lai (2001), where more complete references can be found. 

3.1. Hydrogen Atom 

In a strong magnetic field with 6^1, the electron is confined to the ground Landau 
level ("adiabatic approximation"), and the Coulomb potential can be treated as a 
perturbation. Assuming infinite proton mass (see below), the energy spectrum of the 
H atom is specified by two quantum numbers, (m, z/), where m measures the mean 
transverse separation [Eq. p2|) ] between the electron and the proton, while u specifies 
the number of nodes in the z-wavefunction. There are two distinct types of states in the 
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energy spectrum Em,y The "tightly bound" states have no node in their 2;-wavefunctions 
(z/ = 0). The transverse size of the atom in the (m, 0) state is ^ Pm = [{2m + l)/bY^'^. 
For pm ^ 1, the atom is elongated with L±. We can estimate the longitudinal 

size Lz by minimizing the energy, E ~ — ln{Lz/L±) (where the first term is the 
kinetic energy and the second term is the Coulomb energy), giving 

/ b y 

l^Tr— ^ • (17) 



21n(l/p„) V 2m + 1 
The energy of the tightly bound state is then 



2 



E^o ~ -0.16A l^ln ^—^j (for 2m + 1< 6). (18) 

Here the numerical prefactor comes from solving the Schrodinger equation. The 
coefficient A in ()18p is close to unity for the range of b of interest (1^6^ 10^) 
and varies slowly with b and m (e.g., A ^ 1.01 — 1.3 for m = — 5 when Bu = 1, 
and A ~ 1.02 — 1.04 for m = — 5 when B12 = 10. Note that E^o asymptotically 
approaches —0.5 [ln6/(2m + 1)]^ when b —>■ 00). Obviously the ground-state binding 
energy |i?oo| (e.g., 160 eV at 10^^ G and 540 eV at 10^^ G) is much larger than the 
zero-field value (13.6 eV). For p^ ^ 1, or 2m + 1^6 [but still b ^ (2m + 1)"^ so that 
the adiabatic approximation (|-Emo| ^ b) is valid], we have ~ pH"^, and the energy 
levels are approximated by 

/ b \ 

E^o ~ -0.6 [for 2m + 1 ^ 6 > (2m + 1)-^]. (19) 

\2m +1/ 

Again the prefactor 0.6 comes from solution of the Schrodinger equation. Numerical 
values of E^o for different _B's can be found, for example, in Ruder et al. (1994). Fitting 
formulae for Emo are given in Potekhin (1998) and in Ho et al. (2003). 

Another type of states of the H atom has nodes in the z-wavefunctions {u > 0). 
These states are "weakly bound", and have energies (for b ^ 1) of order E^u ~ — Ry 
(see Lai 2001 and references therein; see also Potekhin 1998 for a comprehensive set of 
fitting formulae). The sizes of the wavefunctions are pm perpendicular to the field and 
Lz ~ z/^ao along the field. 

This simple picture of the H energy levels is modified when a finite proton mass is 
taken into account. Even for a "stationary" H atom, the energy Emu should be replaced 
by Emu + mhucp, where hucp = heB / {nipc) = 6.3 B12 eV is the proton cyclotron energy 
(e.g. Herold et al. 1981). The extra "proton recoil" energy mhucp becomes increasingly 
important with increasing B. Moreover, the effect of center-of-mass motion is non- 
trivial: When the atom moves perpendicular to the magnetic field, an electric field is 
induced in its rest frame and can significantly change the atomic structure, i.e., there 
is a strong coupling between the center-of-mass motion and the "internal" electron 
motion. It is easy to show that even including the Coulomb interaction, the total 
pseudomomentum, 

K = K, + Kp, (20) 
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is a constant of motion. Moreover, all components of K commute with each other. Thus 
it is natural to separate the CM motion from the internal degrees of freedom using K 
as an explicit constant of motion. From Eq. (jZj), we find that the separation between 
the guiding centers of the electron and the proton is directly related to Kj^: 

cB X K 

R-ZC = R-ce — Rep = — • (21) 

For sufficiently small K±, it is convenient to think of the effect of center-of-mass motion 
as the "motional Stark effect", and the kinetic energy associated with the transverse 
motion is Kj_/{2M±rn)- Here the effective mass M±rn depends on the energy levels, and 
increases with increasing B. For large transverse pseudomomentum K^, the structure 
of the "moving" atom is qualitatively different from the stationary atom: The atom 
assumes a decentered configuration, with transverse electron-proton separation ~ Rk, 
and longitudinal separation ~ -R//^, and its energy depends on K± as —R^^ = —b/K_i_ 
rather than the usual dependence (see Lai 2001 and Potekhin 1998 for fitting 
formulae and references). 



3.2. Multielectron Atoms 



The result for H atom can be easily generalized to hydrogenic ions (with one electron 
and nuclear charge Z). The adiabatic approximation (where the electron lies in the 
ground Landau level) holds when p <^ oq/Z, or b ^ Z'^. For a tightly bound state, 
(m, u) = (m, 0), the transverse size is L± ~ pm, while the longitudinal size is 



Zln 



1 



Zpr, 



The energy is given by 



E ~ 



-0.16^^2 




1 2 



(22) 



(23) 



2m + 1^ 

for b ^ (2m + 1)Z^. Results for the weakly bound states {u > 0) can be similarly 
generalized. 

We can imagine constructing a multi-electron atom (with Z electrons) by placing 
electrons at the lowest available energy levels of a hydrogenic ion. The lowest levels 
to be filled are the tightly bound states with u = 0. When oq/Z ^ \/2Z — Ip, i.e., 
b ^ 2Z^, all electrons settle into the tightly bound levels with m = 0, 1, 2, ■ ■ ■ , Z — 1. 
The energy of the atom is approximately given by the sum of all the eigenvalues of 
Eq. ()23j) . Accordingly, we obtain an asymptotic expression for Z ^ \ (Kadomtsev & 
Kudryavtsev 1971) 



E ~ -zHl, 



with Iz = In 



ao 



The size of the atom is given by L± ~ {2Z — 1)^/^ 



P. 



In 



2Z3 



(24) 



ao/{Zl2 



For intermediate-strong fields (but still strong enough to ignore the Landau 
excitation), Z^^^ ^ 6 ^ 2Z'^, many z/ > states of the inner Landau orbitals (states 
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with relatively small m) are populated by the electrons. In this regime a Thomas-Fermi 
type model for the atom is appropriate, i.e., the electrons can be treated as a one- 
dimensional Fermi gas in a more or less spherical atomic cell (Kadomtsev 1970; Mueller 
et al. 1971). The electrons occupy the ground Landau level, with the 2;-momentum 
up to the Fermi momentum pp ^ ^e/h, where is the number density of electrons 
inside the atom (recall that the degeneracy of a Landau level is eB/hc ~ h). The 
kinetic energy of electrons per unit volume is ~ hp^p ~ nl/b"^, and the total kinetic 
energy is Ek ~ R^nl/b'^ ~ Z^/b'^R^, where R is the radius of the atom. The potential 
energy is Ep ~ —Z'^/R. Therefore the total energy of the atom can be written as 
E ~ Z^/{b'^R^) - Z'^/R. Minimizing E with respect to R yields 

Rr^Z'/^-^/', Er.-Z^H^I\ (25) 

For these relations to be valid, the electrons must stay in the ground Landau level; this 
requires Z/R<^ hwce — b, which corresponds to 6 ^ Z^l^. 

Rehable values for the energy of a multi-electron atom for 6 » 1 can be calculated 
using the Hartree-Fock method, which takes into account electron-electron direct 
and exchange interactions in a self-consistent manner. The Hartree-Fock method is 
approximate because electron correlations are neglected. Due to their mutual repulsion, 
any pair of electrons tend to be more distant from each other than the Hartree-Fock wave 
function would indicate. In zero-field, this correlation effect is especially pronounced for 
the spin-singlet states of electrons for which the spatial wave function is symmetrical. In 
strong magnetic fields, the electron spins are in the ground state all aligned antiparallel 
to the magnetic field, the spatial wavefunction is antisymmetric with respect to the 
interchange of two electrons. Thus the error in the Hartree-Fock approach is expected 
to be significantly smaller than the 1% accuracy characteristic of zero-field Hartree-Fock 
calculations (Neuhauser et al. 1987; Schmelcher et al. 1999). 

Accurate energies of He atom as a function of B in the adiabatic approximation 
(valid for h ^ Z"^) were obtained by Virtamo (1976) and Proschel et al. (1982). This 
was extended to Z up to 26 (Fe atom) by Neuhauser et al. (1987) (see also Miller & 
Neuhauser 1991; Mori & Hailey 2002). Numerical results can be found in these papers. 
Neuhauser et al. (1987) gave an approximate fitting formula, E ~ —1&QZ^I^B^2 eV, 
for 0.5 ^ Bi2 ^ 5 (Comparing with the numerical results, the accuracy of the formula is 
about 1% for Z ~ 18 — 26 and becomes 5% for Z ~ 10.) For the He atom, more accurate 
results (which relax the adiabatic approximation) are given in Ruder et al. (1994), Jones 
et al. (1999) (this paper also considers the effect of electron correlation) and in Al-Hujaj 
& Schmelcher (2003a,b). 

Other calculations of heavy atoms in strong magnetic fields include Thomas- Fermi 
type statistical models (see Fushiki et al. 1992; Lieb et al. 1994a,b) and density functional 
theory (Jones 1985, 1986; Kossl et al. 1988; Relovsky & Ruder 1996). The Thomas- 
Fermi type models are useful in establishing asymptotic scaling relations, but are not 
adequate for obtaining accurate binding energy and excitation energies. The density 
functional theory can potentially give results as accurate as the Hartree-Fock method 
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after proper calibration is made (see Vignale & Rasolt 1987,1988). Numerical results of 
He, C, Fe (and the associated ions) for a wide range of field strengths (10^^-10^^ G) are 
given in Medin & Lai (2006). 

The effects of center-of-mass motion on multi-electron systems (heavy atoms and 
molecules) in strong magnetic fields have not been studied numerically, although 
many theoretical issues are discussed in Johnson et al. (1983) and Schmelcher et 
al. (1988,1994). 

While for neutral atoms the center-of-mass motion can be separated from the 
internal relative motion, this cannot be done for ions (Avron et al. 1978). Ions undergo 
collective cyclotron motion which depends on the internal state. However, the existence 
of an approximate constant of motion allows an approximate pseudoseparation up to 
very high fields (see Baye & Vincke 1998). Some numerical results for He"*" moving 
in strong magnetic fields are obtained by Bezchastnov et al. (1998) and Pavlov & 
Bezchastnov (2005). 

3.3. Molecules 

In a strong magnetic field, the mechanism of forming molecules is quite different from 
the zero-field case (see Ruderman 1974; Lai et al. 1992). Consider hydrogen as an 
example. The spin of the electron in a H atom is aligned anti-parallel to the magnetic 
field (flipping the spin would cost huce), and therefore two H atoms in their ground 
states (m = 0) do not bind together according to the exclusion principle. Instead, one 
H atom has to be excited to the m = 1 state. The two H atoms, one in the ground state 
(m = 0), another in the m = 1 state then form the ground state of the H2 molecule by 
covalent bonding. Since the "activation energy" for exciting an electron in the H atom 
from the Landau orbital m to (m -|- 1) is small [see Eq. (fTHj) ]. the resulting H2 molecule 
is stable. Similarly, more atoms can be added to form H3, H4, ■ ■ -. 

The size of the H2 molecule is comparable to that of the H atom. The interatomic 
separation Oeq and the dissociation energy D of the H2 molecule scale approximately as 

aeq-j^, (In 6)2, (26) 

although D is numerically smaller than the ionization energy of the H atom. 

Another mechanism of forming a H2 molecule in strong magnetic fields is to let 
both electrons occupy the same m = Landau orbital, while one of them occupies the 
tightly bound z/ = state and the other the v = 1 weakly bound state. This costs no 
"activation energy" . However, the resulting molecule tends to have a small dissociation 
energy, of order a Rydberg. One can refer to this electronic state of the molecule as 
the weakly bound state, and to the states formed by two electrons in the i/ = orbitals 
as the tightly bound states. As long as In 6 1, the weakly bound states constitute 
excited energy levels of the molecule. 

In the Born-Oppenheimer approximation (see Schmelcher et al. 1988,1994 for 
a discussion on the validity of this approximation in strong magnetic fields), the 
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interatomic potential U{a,R±) is given by the total electronic energy E{a,R±) of 
the system, where a is the proton separation along the magnetic field, and R± is 
the separation perpendicular to the field. Once E{a, R±) is obtained, the electronic 
equilibrium state is determined by locating the minimum of the E{a, R±) surface. [For 
a given a, E{a, R±) is minimal at R± = 0]. The energy curve E{a, R±) can be obtained 
from Hartree-Fock calculations. For large a, configuration interaction must be taken 
into account in the Hartree-Fock scheme (Lai et al. 1992). 

Molecular configurations with R± ^ correspond to excited states of the molecules. 
To obtain E{a,R±), mixing of different m-states in single-electron orbital needs to be 
taken into account. Approximate energy surfaces E(a, R±) for both small R± and large 
R± have been computed by Lai & Salpctcr (1996). 

Numerical results of E{a, 0) (based on the Hartree-Fock method) for both tightly 
bound states and weakly bound states are given in Lai et al. (1992) and Lai & Salpeter 
(1996). Quantum Monte Carlo calculations have also been performed, confirming the 
validity of the method (Ortiz et al. 1995). For example, the dissociation energy of H2 
(neglecting the zero-point energy of the protons) in the ground state is 46 eV for B12 — 1 
and 150 eV for B12 — 10. By contrast, the zero-field dissociation energy of H2 is 4.75 
eV. 

For the ground state of H2, the molecular axis and the magnetic field axis coincide, 
and the two electrons occupy the m = and m = 1 orbitals, i.e., (mi, 7722) = (0,1). 
The molecule can have different types of excitation levels (Lai & Salpeter 1996): (i) 
Electronic excitations. The electrons occupy orbitals other than (mi,m2) = (0,1), 
giving rise to the electronic excitations. The energy difference between the excited 
state (mi, 7712) (with ui = 1^2 = 0) and the ground state (0,1) is of order In 6, as in 
the case for atoms. Another type of electronic excitation is formed by two electrons 
in the (m, u) — (0, 0) and (0, 1) orbitals. The dissociation energy of this weakly 
bound state is of order a Rydberg, and does not depend sensitively on the magnetic 
field strength, (ii) Aligned vibrational excitations. These result from the vibration of 
the protons about the equilibrium separation Ocq along the magnetic field axis, (iii) 
Transverse vibrational excitations. The molecular axis can deviate from the magnetic 
field direction, precessing and vibrating around the magnetic axis. Such an oscillation 
is the high-field analogy of the usual molecular rotation; the difference is that in strong 
magnetic fields, this "rotation" is constrained around the magnetic field fine. Note 
that in a strong magnetic field, the electronic and (aligned and transverse) vibrational 
excitations are all comparable. This is in contrast to the zero-field case, where we have 

Molecules of heavy elements in strong magnetic fields have not been systematically 
investigated until recently. In general, we expect that, as long as a^/Z ^ {2Z — l)^^^p, 
or 6 ^ 2Z^, the electronic properties of the heavy molecule is similar to those of H2. 
When the condition b 3> 2Z^ is not satisfied, the molecule should be quite different and 
may be unbound relative to individual atoms (e.g., at B — 10^^ G bound Fe molecule 
does not exist). Some Hartree-Fock results of diatomic molecules (from H2 up to C2) 
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at 6 = 1000 are given in Demeur et al. (1994). A recent study of molecules of H, He, 
C and Fe using density functional theory for a wide range of field strengths is given in 
Medin & Lai (2006). 

As more atoms are added to a molecule, the energy per atom in a molecule saturates, 
becomes independent of the number of atoms in the molecule. We then essentially have 
a one- dimensional condensed matter (see the next subsection). 



3.4- Condensed Matter 

The binding energy of magnetized condensed matter at zero pressure can be estimated 
using the uniform electron gas model (e.g., Kadomtsev 1970). Consider a Wigner-Seitz 
cell with radius = Z^^^rg (r^ is the mean electron spacing); the mean number density 
of electrons is n^. = Z/(47rrf/3). The electron Fermi momentum pp is obtained from 
rig = [eB /hc){2pp /h). When the Fermi energy p\/{2me) is less than the cyclotron 
energy hujce-, or when the electron number density satisfies 

ne<nB = -i^-r- = 0.0716 (27) 

(or Tj > = 1.49 Z^/^6^^/^), the electrons only occupy the ground Landau level. The 
energy per cell can be written as 

where the first term is the kinetic energy and the second term is the Coulomb energy. 
For a zero-pressure condensed matter, we require dEg/dvi = 0, and the equilibrium Vi 
and energy are then given by 

r,,o ^ 1.90^^/^2/5, (29) 
E,,o ^ -0.395 Z9/562/5_ (3Q^ 

The corresponding zero-pressure condensation density is 

Ps,o ~ 561 A Z-^/^B%^ g cm-3 (31) 

Note that for 6^1, the zero-pressure density is much smaller than the "magnetic" 
density defined in Eq. (j2ZI), i.e., ps,o/ Pb = ('"s/'^i.o)^ = 0.48^^/^6^^/^°. The uniform 
electron gas model can be improved by incorporating the Coulomb exchange energy 
and Thomas- Fermi correction due to nonuniformity of the electron gas (see Lai 2001). 

Although the simple uniform electron gas model and its Thomas-Fermi type 
extensions give a reasonable estimate for the binding energy for the condensed state, 
they are not adequate for determining the cohesive property of the condensed matter. 
The cohesive energy Qg is the difference between the atomic ground-state energy and 
the energy per atom of the condensed matter ground state. One uncertainty concerns 
the lattice structure of the condensed state, since the Madelung energy can be quite 
different from the Wigner-Seitz value [the second term in Eq. (j^H|) ] for a non-cubic 
lattice. In principle, a three-dimensional electronic band structure calculation is needed 
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to solve this problem. So far the only attempt to this problem has been the preliminary 
calculations by Jones (1986) for a few elements and several values of field strengths using 
density functional theory. 

Three-dimensional condensed matter can be formed by placing a pile of parallel 
chains together. The energy difference AEg = |-Es,o| ~ \Eoo\ between the 3d condensed 
matter and the Id chain must be positive and can be estimated by calculating the 
interaction (mainly quadrupole-quadrupole) between the chains. Various considerations 
indicate that the difference is between 0.4% and 1% of \Eao\ (Lai & Salpeter 1997). 
Therefore, for light elements such as hydrogen and helium, the binding of the 3d 
condensed matter results mainly from the covalent bond along the magnetic field axis, 
not from the chain-chain interaction. The binding energies of ID chain for some elements 
have been obtained using Hartree-Fock method (Neuhauser et al. 1987; Lai et al. 1992; 
Lai 2001). Density functional theory has also been used to calculate the structure of 
linear chains in strong magnetic fields (Jones 1985; Relovsky & Ruder 1996; Medin & 
Lai 2006). 

Numerical calculations carried out so far have indicated that for B12 = 1 — 10, linear 
chains are unbound for large atomic numbers Z > 6 (Jones 1986; Neuhauser et al. 1987; 
Medin & Lai 2006). In particular, the Fe chain is unbound relative to the Fe atom; this 
is contrary to what some early calculations (e.g.. Flowers et al. 1977) have indicated. 
Therefore, the chain-chain interaction must play a crucial role in determining whether 
the three dimensional zero-pressure Fe condensed matter is bound or not. The main 
difference between Fe and H is that for the Fe atom at B12 ~ 1, many electrons are 
populated in the 7^ 1 states, whereas for the H atom, as long as 6 ^ 1, the electron 
always settles down in the z/ = tightly bound state. Therefore, the covalent bonding 
mechanism for forming molecules is not effective for Fe at B12 ~ 1. However, for a 
sufficiently large B, when ao/Z ^ ^/2Z + Ip, or B12 ^ 100(Z/26)^, we expect the Fe 
chain to be bound in a manner similar to the H chain or He chain (Medin & Lai 2006). 
The cohesive property of magnetized condensed matter is important for understanding 
the physical condition of the "polar gap" of the neutron star (see §9). 

4. Propagation of Photons in a Magnetized Plasma 

The magnetized plasma around a neutron star is anisotropic and birefringent, and 
significantly influences the polarization state of a photon. Many of the radiative 
processes discussed in ^ are polarization-dependent (e.g., photons in different 
polarization modes often have very different opacities). Thus it is important to 
understand the propagation of photons in a magnetized plasma. The atmosphere of a 
neutron star (with density 10~^ — 10^ g/cm^) can be characterized as a cold plasma for 
most purposes, while the much more tenuous magnetosphere likely comprises relativistic 
electron-positron pairs. Even in a pure vacuum, the effect of quantum electrodynamics 
induces birefringence through vacuum polarization. 
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4.1. Dielectric Tensor of Cold Plasma 

Following Ginzburg (1970), we consider a cold plasma composed of electrons and ions 
(with charge, mass and number density given by — e, mg, Ug and Ze, mi = Anip, Ui = 
n^/Z, respectively; here Z is the charge number and A is the mass number of the ion) in 
an external magnetic field B. The electrons and ions are coupled by collisions, with the 
collision frequency u^i. We can generalize Ginzburg (1970) by including the radiative 
dampings of electrons and ions, with the damping frequencies Ure and I'ri, respectively. In 
the presence of an electromagnetic wave with the electric field E oc e~*'^*, the equations 
of motion for a coupled electron-ion pair are: 

" (32) 



meFe = — eE — 

c 
Ze 

miici = ZeEi H rj x B 

c 



X B - mel/ej(re - i-j) - meVre^e, 



m^Zu^Aii 



(33) 



Note that the drag coefficient due to e-ion collision for electron is me{(JeiVei)ni = mgUei 
(where is the cross-section, Vei is the relative velocity), while the drag coefficient 
for ion is me{<JeiVei)ne = meZvei since Ue = Zui. Solving Eqs. (jS2I)-(|Sni) yields 
r = Te — Tj, and the polarization of the plasma P = UiZeri — ngere = — neCr. The 
electric displacement vector is then D = E -|- 47rP = e*^^^ ■ E. In the coordinate system 
XYZ with B along Z, the plasma dielectric tensor is given by 
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In Eqs. (IHHjl and (jHEl), we have defined the dimensionless quantities 
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where Uce = eB/{mec) is the electron cyclotron frequency, ujd = ZeB/{mic) is the 
ion cyclotron frequency, Upe = {'inriee^ /meY/"^ is the electron plasma frequency, and 

lei = i^ei/^, Ire = ^re/^, and 7^ = ^ri/^ are given by 



[A-KUiZ'^e^ /miY^"^ is the ion plasma frequency. The dimensionless damping rates 
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where g^ is the Gaunt factor (see Potekhin & Chabrier 2003 and references therein). 
Including these damping terms in the dielectric tensors allows one to obtain the 
appropriate expressions for the radiative opacities using the imaginary parts e^^\ If 
we neglect damping {■jre = Iri = Jei = 0), Eqs. (jHSj) and (jHUj) reduce to 



1/2 



1/2 



1 - 1 



Ui 



1 — Up 1 — Uj 



1] = 1 - Ve - Vi. (41) 



4-2. Vacuum Polarization 

In strong magnetic fields, vacuum polarization can significantly influence photon 
propagation. According to quantum electrodynamics, a photon may temporarily 
convert into virtual electron-positron pairs, and since the pairs are "polarized" by 
the external magnetic field, the magnetized vacuum has a nontrivial dielectric tensor 
e and a nontrivial permeability tensor fj, (e.g., Heisenberg & Euler 1936, Schwinger 
1951; Adler 1971; Tsai & Erber 1975; see Schubert 2000 for extensive bibliography). 
In the low-energy limit, huj <^ meC^, the Euler-Heisenberg effective Lagrangian can 
be applied. The relevant dimensionless magnetic-field parameter is (3 = B/Bq, where 
Bq = mlc^/eh = 4.414x10^^ G. The vacuum dielectric tensor e and inverse permeability 
tensor fj,^ = fx take the form e = I + Ae^^^ and p, = I + Afj,^""^ (where I is the unit 
tensor), with 

Ae(^) =al + qBB, (42) 
A/i^'^) = al + mBB, (43) 

where a, q, and m are functions of (3. In the limit /3 -C 1, the vacuum polarization 
coefficients are given by (Adler 1971) 
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where a = /he = 1/137 is the fine-structure constant. For arbitrary B, the vacuum 
polarization coefficients have been obtained by Heyl & Hernquist (1997) in terms of 
special functions and by Kohri & Yamada (2002) numerically. Convenient expressions 
are given by Potekhin et al. (2004): 
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where X(^) (with ^ = P ^) is expressed through the Gamma function T{x): 
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The following simple fitting formulae can be used for all values of /3: 
2a / (3^ 1 + 0.25487 

-- (50) 



^ 4571^^ 1 + 1.33^ + 0.56/32' 

"^"^ ~3^ 3.75 + 2.7/35/4 + /52- (^^) 
Equations ()49|) - (jHT|) exactly recover the weak-field limits ()44j) and the leading terms in 
the high-field {b ^ 1) expansions (Eqs. [2.15]-[2.17] of Ho & Lai (2003). The maximum 
errors are 1.1% at /3 = 0.07 for 2.3% at /3 = 0.4 for and 4.2% at /5 = 0.3 for 
(j^Tj) . The derivatives of a, q, m with respect to (3 have similar accuracies. 

^^.5". Photon Modes in Cold Plasma with Vacuum Polarization 

Including vacuum polarization, the dielectric tensor e and inverse permeability tensor 
fi~^ = fj, can be written as e = e*^^) + Ae*^''^ = I + Ajl^'"'^ (where I is the unit tensor). 
Thus the dielectric tensor of the combined plasma+vacuum medium still be written in 
the form of Eq. except e^e' = e + a, r]^r]' = r] + a + q. 

An electromagnetic (EM) wave with a given frequency uj satisfies the wave equation 

V X (/X • V X E) = — e ■ E. (52) 

For normal modes propagating along the 2;-axis, with E oc e*'^='=^, this reduces to 

- klz x[fi-{zx E±)] = — e ■ E±, (53) 

where the subscripts "±" specify the two modes ("plus- mode" and "minus-mode"). 
In the xyz coordinates with k along the ^-axis and B in the x-z plane (such that 
B X k = sin 9b y, where 9b is the angle between k and B), we write the electric field of 
the mode as Ei± oc {iK±, l,iKz±), where the ellipticity K± = —lE^/Ey is given by 

K^ = (3,±pl + r, (54) 

with r = 1 + {m / a) SIT? 9 B — 1 (where a = 1 + d) and the (complex) polarization 
parameter /3p is 

e'"^ - g"^ - e't]' {I + m/a)siv? 9b .r■r■^ 
2gri' cos 9 b 

The refractive index n± = ck±/uj is given by 

n| = ^f- + ^cos^BV (56) 
ae33 ^± / 

where £33 = e' siv? 9b + rj' cos^ 9b- 

The polarization parameter j3p directly determines the characteristics of photon 

normal modes in the medium. For fg ^ 1? the real part of (3p can be written as 

Re(/3p) = [3q(3v, where 

ul'"^ sir? 9b ,^ . 
2 cos9b 
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-.=(-^rW) > ^ii = ( . I • (60) 



and 

Pv^l+ (^ + ^)(^--^) . (58) 

UeVe 

In general, the modes are elliptically polarized, and the ellipticity depends on Pp. 

A special case is the magnetized vacuum with zero plasma density. The 
extraordinary mode (also called ±-mode) and 0-mode (also called ||-mode) are simply 

E^ = (0,1,0), E||oc(l,0,ir.) with K^^ I^^^^BCOse^ ^ ^^g^ 

a + gcos^^t^B 

The indices of refraction are 

1/2 ^ ^ / a + q V^' 
^a + msin^OB/ ' " \a + gcos^^B, 
Clearly, acting by itself, the birefringence from vacuum polarization is significant (with 
the index of refraction differing from unity by more than 10%) for B > 300Bq — for 
such superstrong field strengths, one can expect appreciable magnetic lensing effect from 
photons emitted near the neutron star surface (e.g., Shaviv, Heyl & Lithwick 1999). 

Now consider the case of finite plasma density. For concreteness, let us consider 
the regime ^ 1, i.e., the photon energy e <^ E^e- This is relevant for propagation 
of thermal photons in neutron star atmospheres. Clearly, for generic values e and 
6b, 1/5 1 ^ 1, the two modes are almost linearly polarized: the extraordinary mode 
has \K\ -C 1, and its E is mostly perpendicular to the fc-B plane; the ordinary 
mode has li^l ^ 1, and is polarized along the A;-B plane. This distinction of the 
two modes manifests most significantly when we consider how they interact with 
matter: the ordinary-mode opacity is largely unaffected by the magnetic field, while 
the extraordinary-mode opacity is significantly reduced (by a factor of order uj'^/uol^ 
from the zero-field value. However, this distinction becomes ambiguous when < 1 
{\K\ ~ 1). Obviously, /5 = for u; = But even for general energies (e ^ hud), a 
photon traveling in an inhomogeneous medium encounters /3 = when the condition 
f e = g + m is satisfied. This is the vacuum resonance (see Gnedin et al. 1978). Using 
htUpe = fi^A-KneC^ / meY^"^ = 28.71 (Y'ePi)^''^ eV, where Y^, is the electron fraction of the 
gas, and pi is the density in unity of 1 g cm^^, we find that for a given photon energy, 
vacuum resonance occurs at the density (Lai & Ho 2002) 

pv = 0.964y-i(5i4ei)V~'g cm-3, (61) 
where ei = e/(l keV) and 

is a slow-varying function of /9 [/ = 1 for /? <^ 1 and / — > (/?/5)^/^ = 0.673 5^4^ for 
/3 ^ 1; / varies from 0.99 at B14 = 1 to 6.7 at B^ = 100]. Qualitatively, the meaning 
of the vacuum resonance is the following (see Fig. 12)): For p > pv, the plasma effect 
dominates the dielectric tensor, while for p < pv, vacuum polarization dominates. Away 
from the resonance, the photon modes (for e <^ Ece) are almost linearly polarized as 
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discussed above. Near p = pv, however, the normal modes become circularly polarized 
as a result of the "cancellation" of the plasma and vacuum effects — both effects tend 
to make the mode linearly polarized, but in mutually orthogonal directions. 

Including the damping terms in the dielectric tensor gives rise to the phenomenon 
of mode collapse (see Soffel et al. 1983; Lai & Ho 2003a). This occurs when the two 
polarization modes become identical, i.e., = K_ or = n_. Obviously, at the 
point of mode collapse, the modal description of radiative transfer breaks down and a 
rigorous treatment of radiative transfer requires solving the transport equations for the 
four photon intensity matrix, or Stokes parameters (see Lai & Ho 2003a). 



4.4- Wave Propagation in Inhomogeneous Medium 

Because of the strong gravity, a neutron star atmosphere has a rather large density 
gradient, with density scale height of order a few em's for T = 10® K. When a 
photon propagates in such an inhomogeneous medium, its polarization state will evolve 
adiabatically (i.e. following the or K_ curve in Fig. 12)) if the density variation is 
sufficiently gentle and the modes are sufficiently distinct (i.e., ni — n2 is sufficiently 
large). Away from the vacuum resonance, the adiabatic condition is easily satisfied. At 
the vacuum resonance, adiabaticity requires (Lai & Ho 2002) 

e ^ ead(5, Ob, Hp) = 1.49 (/ tan^^^ll - u,\f'' ( ^ j keV, (63) 



Hp 



where Hp = \dz/d\np\ is the density scale height (evaluated at p = pv) along the ray. 
For an ionized Hydrogen atmosphere. Hp ~ 2kT / {rupg cos 9) = 1.65 T6/((7i4 cos ^) cm, 
where T = 10® Tg K is the temperature, g = lO^^gi^ cm s~^ is the gravitational 
acceleration, and 6 is the angle between the ray and the surface normal. In general, the 
probability for a nonadiabatic "jump" is 

Pjump = exp [-(7r/2)(e/ead)'] • (64) 

Thus for e ^ 1.3ead, the polarization evolution is highly adiabatic (with Pjump ^ 0.03) 
even at the vacuum resonance. In this case, an X-mode (0-mode) photon will be 
converted into a 0-mode (X-mode) as it traverses the vacuum resonance. This resonant 
mode conversion is analogous to the MSW effect of neutrino oscillation (e.g., Haxton 
1995). One could also say that in the adiabatic limit, the photon will remain in the 
same plus or minus branch, but the character of the mode is changed across the vacuum 
resonance. Indeed, in the literature on radio wave propagation in plasmas (e.g., Budden 
1961; Zheleznyakov et al. 1983), the nonadiabatic case, in which the photon state jumps 
across the continuous curves, is referred to as "linear mode coupling". It is important 
to note that the "mode conversion" effect discussed here is not a matter of semantics. 
The key point is that in the adiabatic limit, the photon polarization ellipse changes 
its orientation across the vacuum resonance, and therefore the photon opacity changes 
significantly (see Lai & Ho 2003a). 
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Figure 2. The polarization ellipticity of the photon mode as a function of plasma 
density near the vacuum resonance. The two curves correspond to the two different 
modes. In this example, the parameters are B — 10^^ G, e — 5keV, Ye = 1, and 
(^kB — 45°. The ellipticity of a mode is specified by the ratio K — —iE^/Ey, where 
Ex (Ey) is the photon's electric field component along (perpendicular to) the k-B 
plane. At densities away from the resonance density py, the two modes are almost 
linearly polarized (with polarization ellipses orthogonal to each other): The 0-mode 
is characterized by \K\ 1, and the X-mode |i^| <C 1. At p = pv, both modes are 
circularly polarized. In the adiabatic limit, an 0-mode (X-mode) photon from the 
high-density region will convert to the X-mode (0-mode) as it traverses the vacuum 
resonance density pv, with its polarization ellipse rotated by 90°. 
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The concept of adiabatic evolution of photon polarization is also important when 
we consider wave/photon propagation in neutron star magnetospheres. For example, 
the normal modes of X-ray photons in the magnetosphere are determined by vacuum 
polarization (since p -C pv)- As a photon (of a given mode) propagates from the polar 
cap through the magnetosphere, its polarization state evolves adiabatically following the 
varying magnetic field it experiences, up to the "polarization-limiting radius" r^i, beyond 
which the polarization state is frozen. This polarization-limiting radius is located where 
the adiabatic condition breaks down (see Heyl et al. 2003; Lai & Ho 2003b). Another 
example concerns the radio wave propagation in the magnetosphere, whose dielectric 
property is dominated by highly relativistic pair plasmas. Such a propagation effect 
may be responsible for some of the observed polarization changes in radio emission of 
pulsars (e.g., Cheng & Ruderman 1979; Barnard 1986; Melrose & Luo 2004). 

4-5. Warm Plasma and Relativistic Plasma 

While the dielectric tensor for cold plasma was derived classically, the quantum 
calculation, incorporating the quantized nature of electron motion transverse to the 
magnetic field, yields the same result when the photon wavenumber /c — >■ (e.g., 
Canuto & Ventura 1972; Pavlov et al. 1980). More precisely, since the size of the 
Landau wavefunction is of order p = {hc/eBY^^ [see Eq. this requires k±p <^ 1 
(where k± is the wave number perpendicular to B, or the photon energy must satisfy 

The cold plasma approximation is valid when the thermal velocity of the electron, 
Vt = {2kBT /meY/"^ ^ is much less than the phase velocity of the wave [uj/k ~ c), and 
oj is not too close to the cyclotron frequency, i.e., for \uj — ujce\ ^ {uj/c)Vt (we neglect 
the ion component here). The cold plasma approximation also neglects electron recoil. 
Taking into account the electron motion, the cyclotron resonance denominator uj — ujce 
in the cold plasma dielectric should be replaced by u; — k\\V\\ — ojce (where k\\ and V\\ are 
the wavenumber and electron velocity parallel to the magnetic field), and the electron 
recoil can also be included by ^ ± Tiky (2me). Then the dielectric tensor should be 
averaged over the electron velocity distribution function. In the very strong magnetic 
field limit where electrons occupy only the lowest Landau level {ksT <^ -Ece), the cold 
plasma expressions PT|) are replaced by (Kirk 1980) 




(65) 



ce 



(66) 
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2VTk\ 
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[the hk±/{2meVT) term arises from electron recoil], and W{z) is the plasma dispersion 
function 



1 f°° e dx , , 



which should be evaluated for Uce i^ce/(l + ^7re) with ■jre > [see Eq. fjHK|) ]. 
Expressions (jUBjl - (jU7jl are valid in the nonrelativistic limit, hu <C mgC^ and ksT <C meC^. 
One can check that for |ci;— ti^ce — ^A;|/(2me)| 3> iVrfcyl, the cold plasma limit is recovered. 
The net effect of finite temperature is that the cyclotron resonance is softened from the 
sharp peak of the cold plasma theory to a broadened one with width Au ~ k\\VT- Also, 
the electron recoil shifts the resonance by Ae {EccCosOb)'^ /{'imec'^), which becomes 
appreciable as B approaches Bq. Thus a quantitative description of the cyclotron 
resonance for B > Bq requires relativistic theory (see 

Similar qualitative behavior appears in the relativistic theory (e.g., Svetozarova 
& Tsytovich (1962); Melrose 1974; Pavlov et al. 1980; Nagel 1981; Meszaros & Nagel 
1985; see Meszaros 1992 for a review). The polarization tensor is directly proportional 
to the forward-scattering amplitude of a photon in the plasma. Quantum states and 
momentum distribution of electrons can be summed or averaged. The advantage of such 
an approach is that the photon polarization modes are calculated in a self-consistent 
manner as the differential scattering cross section, which depends on photon modes (see 

The magnetosphere of a radio pulsar is believed to be filled with ultra-relativistic 
electron-positron pairs (see In the polar-cap region, outwardly streaming pair 

plasmas are produced in an electromagnetic cascade (Sturrock 1971), with bulk Lorentz 
factor 7 > 10^ (e.g., Daugherty & Harding 1982; Zhang & Harding 2000; Hibschman 
& Arons 2001). While the polarization state of high-energy emission is determined by 
the vacuum polarization effect, radio emission generated in the inner magnetosphere 
can be affected when propagating through such a plasma (e.g., Cheng & Ruderman 
1979; Barnard 1986). For example, the observed characteristics of circular polarization 
(Radhakrishnan & Rankin 1990; Han et al. 1998) are thought to develop as a 
propagation effect (e.g., Petrova & Lyubarskii 2000; Melrose & Luo 2004; Petrova 
2006). The most important region is associated with the cyclotron resonance, where 
7a;(l— /3 cos 9b) = uj^e- Dispersion in a pulsar plasma has been studied extensively, under 
various assumptions about the relative electron-positron densities and their momentum 
distribution, which are uncertain (e.g., Arons & Barnard 1986; Lyutikov 1998; Melrose 
et al. 1999; Asseo & Riazuelo 2000). A simple model is based on applying a Lorentz 
transformation to the magneto-ionic (cold plasma) theory (Melrose & Stoneham 1977; 
Melrose & Luo 2004), which ignores the random motions of the particles. Such a model 
has been used to interpret the circular polarization in pulsar radio emission (Melrose & 
Luo 2004). The effect of vacuum resonance in the magnetosphere plasma of pulsars and 
magnetars is studied by Wang & Lai (2006). 
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5. Radiative Processes for Free Electrons 

At low magnetic field strengths, many of the radiative processes we discuss below, 
including cyclotron absorption and emission and Compton scattering, may be described 
by classical physics. In neutron star magnetic fields approaching the critical field, a 
relativistic quantum description is required to accurately compute the rates and the 
photon spectrum of radiation, even for non-relativistic electrons. For these processes, 
we will first review the classical descriptions before discussing the quantum description 
and discuss the regimes in which each is appropriate. Other processes we will discuss, 
including one-photon pair production and annihilation, photon splitting and bound 
pair creation, take place only in strong fields since they do not conserve energy and 
momentum in free space. In a magnetic field, perpendicular momentum of particles 
is not conserved in transitions between Landau states since the field effectively is able 
to supply or absorb momentum. The momentum conservation along the field couples 
directly to the translational invariance of the system (see Section 2) parallel to B. The 
conservation equations for energy and parallel momentum for transitions of electrons or 
positrons between initial state {n,s,p) and final state {n',s',p'), where s = (cTz + l)/2 
is the electron spin state, p = cpz and n labels the Landau level [see Eq. (fTH|l ]. read 

En - = ±e (70) 

p' = PT ecos6' (71) 

resulting in emission (upper sign) or absorption (lower sign) of photons at angle 9 to 
the magnetic field with energy e. [In this Section, in contrast to previous Sections, we 
will express all photon and particle energies in units of mgC^.] The photon energy, as 
determined from the above kinematic equations, is 

{En - p COS 6) - [{En -pcosOf- 2B' [u - n') sin^ 6^/'^ 

^ — . 2 o ' ^ ) 

sm d 

where in this section, we use B' = B/Bq as the dimensionless magnetic field parameter 
(equivalent to the (3 parameter in section H}. Unless the photon is emitted or absorbed 
with angle 6 = 90°, the electron will experience a recoil along the field direction, given 
by Eqn (f7T|) . which is needed to determine its final energy. The electron wavefunctions 
that have been adopted in the literature for studying the processes presented in the 
following subsections were discussed in §21 

5.1. Cyclotron Absorption 

Cyclotron absorption, the inverse of cyclotron emission (see ^5.2p . is a first-order process 
in which a photon excites a particle to a higher Landau state. The classical, non- 
relativistic cross section for absorption of a photon by an electron at rest in the ground 
state (n=0) (or in the rest frame of a moving electron) was derived by Blandford & 
Scharlemann (1976) 

^at(IU) = 2avr^A^5(e-iVS') (^E'sin^^)'' ' ^I^' ^^3) 
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where is the harmonic number, a = 1/137, A = h/{mf,c), 9 is the angle of photon 
propagation to the magnetic field and the quantity in parentheses refers to the two 
incident photon linear polarizations (||, _L) or (0,X) modes, with electric vectors parallel 
and perpendicular to the plane containing the photon wavevector and the magnetic field. 
The 5-function restricts the energy of the absorbed photon to be a harmonic, = NB', 
of the cyclotron energy = Ece/meC^ = B' . 

The relativistic QED cyclotron absorption cross section was first derived by 
Daugherty & Ventura (1978) for electrons initially in the ground state. The required 
energy for excitation to state now follows from the relativistic kinematic equations 
dZOl) and dZH), setting E„ = 1, n = 0, = iV and p = 0, 

= [(1 + 2NB' sin' Of'^ - 1]/ sin' 6 (74) 

for a photon propagating at angle 9 to the field. Because the electron experiences a 
recoil of p' = eNCOs9 on absorption, the cyclotron harmonics are actually anharmonic 
in high magnetic fields, so that the energy difference between successive harmonics will 
decrease. The cyclotron absorption rate, summed over final spin states of the electron, 
is (Harding & Daugherty 1991) 



where 



En ^(iV-1)! 
£^ sin^ 9 



Z 

[cos'^9,1) + ^(sin^^, 0) 



(75) 



(76) 



2B' ' 

The ^-function requires that the absorbed photon energy have only the value given by 
equation and effectively makes only a function of the incident photon angle 9. 
The above expression reduces to the classical limit of equation (fTSjl when NB' << 1, 
where e ~ NB', Z ^ N'^B' sir? 9/2 <C 1 and ~ 1. Therefore, the relativistic formula 
(|75|) should be used when NB' ~ 0.1, that is, when electron recoil becomes significant, 
which is the case in fields above ~ 10^^ Gauss. 



5.2. Cyclotron and Synchrotron Radiation 

Cyclotron radiation is the inverse of the cyclotron absorption process discussed above, 
and results from downward transitions between Landau levels. The classical cyclotron 
radiation formula was first presented by Schott (1912) and assumes a continuous circular 
orbit of the particle in a magnetic field with negligible energy loss. The emission is 
characterized by the particle velocity /3 and pitch angle, ip = ta.'n~^{P±/ where P\\ 
and P± are the velocity parallel and perpendicular to the field direction. If the particle 
energy is above the cyclotron energy, the radiation is the sum over a number of harmonics 
of the fundamental cyclotron energy. The radiated power at frequency u and angle 9 to 
the magnetic field is given by (Sokolov & Ternov 1986) 

,y . /' „ T"J i'^^)' m)^AJl\()\(77) 

c ^ JO (1 - Ml cost^)'^ I V smfc' / I 
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where cu — eB/{'ymec) is the particle gyrofrequency and 

uP±_ sin 9 

^ = l-/3||Cose- ^^^^ 
The polarized form of equation (77), in the case where j3\\ = 0, is given in Sokolov and 
Ternov (1986). The spectrum, given by Bekefi (1966), Canuto and Ventura (1977) and 
Brainerd and Lamb (1987), display characteristic harmonics at low energy. When the 
particle energy is relativistic and the emission is dominated by high harmonics, it is 
called synchrotron radiation. The formula becomes somewhat simpler and the emission 
power per unit photon energy is (e.g. Jackson 1975) 

where cub = <^ce = eB/{meC) is the cyclotron frequency, T = 7_B', 7 is the electron 
energy and the k function is 

poo 

<y) = y / K,/s{x)dx. (80) 

The characteristic energy of a synchrotron photon is e ~ 7^5' sin ^. 

In high fields approaching the critical field, the classical description should be 
replaced with a relativistic quantum description to accurately compute the radiative 
rates and photon spectrum, even for non-relativistic electrons. The classical value of 
the critical radiation frequency, j^eB, exceeds the electron kinetic energy, (7 — 1) when 

^'er^ > 1- (") 

The classical emissivity therefore violates conservation of energy. The classical formula 
also overestimates the spectral emissivity when 

Tsin^ = > 0.1. (82) 

There is a consequent reduction of the energy loss rate for large values of T (Erber 
1966). Finally, taking into account electron recoil gives harmonic energies which are 
not simple multiples of the cyclotron frequency, as in the non-relativistic case. The 
polarization and spin dependent transition rates for quantum synchrotron radiation are 
given by Sokolov & Ternov (1968, see also Harding & Preece 1987) and should be used 
for transitions between low Landau states in field B' ~ 0.1. The exact QED spectrum 
for transition from state n to the ground state n' — 0, averaged over electron spin, is 
(Latal 1986, Baring et al. 2005) 



de A B'€lT{n + l) \e. 



- 1 



{en - e)nB' - e ( ^\ /qqn 

: exp —n — (83) 



^(£„-l-e)(£„ + l-e) ^ V 

where e„ = ^/l + 2nB' and e_ = uB' /sn- The rate is more complicated for transitions 
between higher Landau states, involving Laguerre polynomials. The cyclotron decay rate 
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Figure 3. Spin-averaged cyclotron decay rates (solid lines), in units of the cyclotron 
frequency lob, from Landau state n as a function of magnetic field strength B in units 
of the critical field Bq . Also plotted are the corresponding non-relativistic decay rates 
(dashed hues). 



from state n is found by integrating over the radiated photon energy in each transition, 
Tnn' and then summing over the final states n': 

n-l 

Tn = E Tnn'. (84) 
n'=0 

The spin-averaged decay rates for some selected states, as a function of field strength, 
are shown in Figure El and compared to the classical (non-relativistic) values (Herold 
et al. 1982). One can see that the QED decay rate departs from the classical rate 
for B'n ~ 0.1. The spin-dependent decay rates have been discussed by Daugherty & 
Ventura (1978) and by Melrose & Zheleznyakov (1981) for the non-relativistic case, and 
by Herold et al. (1982) for the relativistic case. Generally, the probability for spin-flip 
transitions, those in which the electron changes its spin state, is lower than for non 
spin-flip transitions. The ratio of the spin-flip decay rate to the non-spin-flip decay rate 
is oc B'^^ for the n = 1 — transition. 

An asymptotic form of the spectrum, averaged over spin and polarization, valid 
for relativistic elections, transitions between high Landau states (n, n' ^ 1) and zero 
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longitudinal momentum {ip = 90°), is (Sokolov & Ternov 1968) 



W = (1 - 1) .(,) + y^hni - i)K,M 



(85) 




(86) 



37T[l-(e/7)] 



One can see that the above expression incorporates a kinematic cutoff in the spectrum 
at e = 7 that avoids energy violation. When e <C 7 and 5' ^ 1, formula (jHH|) reduces 
to the classical formula (f?^ . 

There are several notable features of cyclotron/synchrotron radiation in high 
magnetic fields. Classically and at low field strengths, rates are highest for single 
harmonic number transitions and decrease monotonically with increasing harmonic 
number, with transitions from high Landau states to the ground state being very 
improbable. But for field strengths B' ~ 0.2, transitions from high to low states become 
more probable and transitions to the ground state can actually dominate over even 
single harmonic number transitions (White 1974, Harding Sz Preece 1987). The result 
is that electrons in very high magnetic fields radiate energy not in small steps but often 
in one large transition, emitting a photon equal to its kinetic energy. The spectrum 
displays an enhancement just before the kinematic cutoff, as shown in Figure |3 

The energy loss rate of an electron emitting synchrotron radiation, resulting from 
integration of the spectrum of Eqn ()85|) over photon energy, is shown in Figure El as 
a function of the parameter T. For low values of T <^ 0.1, the synchrotron loss rate 
follows the well-known classical dependence 7 oc T^, but begins to depart from the 
classical dependence for values of T ~ 0.1. At large values of T ^ 1, the loss rate 
follows a milder dependence 7 oc T^/^. This softening of the loss rate occurs because 
the total rate of emission declines at high T (Baring 1988). 

5.3. Compton Scattering 

In neutron star magnetic fields, the cyclotron decay rate is high enough so that 
nearly all particles occupy the ground state. The cyclotron de-excitation rate rcyc ~ 
3 X 10^^Bf2 s~^, is much larger the collision rate, rcoii ~ 5 x 10^(ne/10^^ cm~^) Bi2^^ 
(Bonazzola et al. 1979). Radiative (rather than collisional) processes thus control 
the Landau state populations. Radiative transitions thus dominate over collisions in 
astrophysical sources having strong magnetic fields. An electron that absorbs a cyclotron 
photon will almost always de-excite by emitting another photon, rather than collisionally 
de-exciting, so that the end result of the process is a scattering of the photon rather than 
a true absorption. This means that in strong fields resonant scattering dominates over 
absorption. Cyclotron absorption is most accurately treated as a second-order process, 
as resonances in the Compton scattering cross section, where the excited state of the 
particle is a virtual state for which energy and momentum are not strictly conserved. 
The cyclotron lines are thus broadened by the intrinsic width which is equal to the 
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Figure 4. Quantum synchrotron spectra for a particle with Lorentz factor 7 = 100, 
initial pitch angle, sinip = 1 and different field strengths B' in units of the critical 
field, computed from the asymptotic formula in Eqn 

inverse of the decay rate from that state (Harding & Daugherty 1991, Graziani 1993), 
although Doppler broadening usually dominates at most angles. The addition of the 
decay rate to the resonant denominator also renders the cross section finite at resonance. 
The cross section for Compton scattering in a magnetic field was first studied in the 
non-relativistic limit by Canuto, Lodenquai & Ruderman (1971) and the full QED cross 
section has been computed by Herold (1979), Daugherty & Harding (1986) and Bussard, 
Meszaros & Alexander (1986). Since the non-relativistic treatment is limited to dipole 
radiation, only scattering at the cyclotron fundamental is allowed. In the relativistic 
(QED) treatment scattering at higher harmonics is allowed, including Raman scattering, 
in which the state of the particle after scattering is higher than the initial state. 

In free space, the total electron scattering cross section is just the Thomson cross 
section, ctt = 6.6 x 10~^^cm^, for photon energies which are non-relativistic in the 
electron rest frame, or when €7 -C m^. At higher energies, relativistic effects are 
important both in the kinematics and in the cross section (e.g. Rybicki and Lightman 
1979). The photon energy change in the electron rest frame, due to electron recoil, can 
no longer be ignored and the Klein- Nishina cross section (e.g. Jauch & Rohrlich 1980) 
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is appropriate in the relativistic (QED) regime. 

The classical, non-relativistic (Thomson) limit of the magnetized scattering cross 
section (Canuto et al. 1971, Ventura 1979) has a strong dependence on photon frequency, 
angle to the magnetic field and polarization. The non-relativistic total scattering cross 
section in the electron rest frame for linearly polarized photons takes the form (Blandford 
and Scharlemann 1976): 



(T|| = (Jt { sin^ ^ 2 ^ 



e2 



+ 



tT± = ^ r ' + ^ I . (87) 

for polarization states (|| , -L) where to the photon electric vector parallel or perpendicular 
to the plane formed by the photon wavevector and the field. Here 6 and e are the angle 
and energy of the incident photon with respect to the field in the electron rest frame, 
and es = TieB / {mec) is the cyclotron energy. 

The main effects of the magnetic field on electron scattering is the appearance 
of the cyclotron resonance and a strong dependence of the cross section on photon 
energy and incident angle. For photon energies well above the resonance, e ^ es, 
a = aT for both polarizations. For photon energy below the resonance, e -C e^, 
o"|| ~ crr[sin^6' + cos^ 9{e/eBY] and a± ~ o"r(e/eB)^. A magnetized plasma becomes 
quite optically thin for propagation parallel to the magnetic field for photon frequencies 
below the cyclotron frequency. 
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As in the case of cyclotron absorption, the non-relativistic cross section for 
scattering is not accurate for B' ~ 0.1 and furthermore describes scattering only in 
the fundamental and so cannot be used to treat scattering in higher harmonics. As we 
will discuss in more detail in ^IH cyclotron line modeling at low field strengths is still 
possible using the absorption cross section for harmonics above the fundamental, where 
the scattering is treated as cyclotron absorption followed by cyclotron emission of the 
photon. Even at high field strengths, use of the relativistic absorption cross section can 
give a reasonable approximation to the scattering cross section for low harmonics. But 
this treatment neglects the contribution from non-resonant scattering which is important 
for photon energies away from resonance. 

The relativistic magnetic scattering cross section is quite complicated even for the 
simplest case of ground state-to-ground state scattering (Herold 1979), and is even 
more unwieldy for the case of scattering to arbitrary Landau states (Daugherty Sz 
Harding (1986) and Bussard, Meszaros & Alexander 1986). However, a number of 
interesting features appear in the relativistic scattering cross section (Daugherty & 
Harding 1986). One is the possibility of ground state-to-ground state scattering through 
higher intermediate Landau states (i.e. n 0, producing resonant contributions 
to the cross section at higher harmonics. Another possibility is a — — n scattering 
where a photon of energy e ~ nB' is almost completely 'absorbed' producing a scattered 
photon having very low energy. The inverse process in which an electron in an excited 
state scatters to the ground state, can convert an X-ray photon to a 7-ray photon 
(Brainerd 1989). Another possibility that arises solely in QED is two-photon emission, 
where an electron makes a transition between two Landau states with the emission of 
two photons (Alexander & Meszaros 1991a, Semionova & Leahy 1999, Melrose & Kirk 
1986) Gonthier et al. (2000) have derived a simplified analytic approximation to the 
relativistic Compton scattering cross section for the case where the incident photon is 
parallel to the field, in which case there is a resonance only at the fundamental and 
scattering to higher Landau states can effectively be neglected. Such a situation would 
apply to scattering by a relativistic electron moving along the magnetic field, where 
photons would primarily appear predominatly in a Lorentz cone (I/7) centered on the 
field direction in the electron rest frame. The approximate expression for the total 
relativistic, polarization-dependent scattering cross section for the case 6* = 0° is 




3(Tt 



{f(e)-2eh{e)} 
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where 




+ e(l + 2e) ln(l + 2e) + 2e. 
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/(e) = - ln(l + 4e) + e(l + 2e) ln(l + 2e), and 
h{e) = 



(89) 



tan 



1+e 



In 



(l+e+^e(e-2)y 



l+4e 



for e < 2, and 
for e > 2. 



This expression reduces to the nonrelativistic hmit for small e. Figure El shows this 
approximation compares to the exact QED cross section for the same case of incident 
photon angle sin 6' = 0. For incident photon energies above the cyclotron energy, the 
electron can be excited to higher Landau states and this contribution (for states up 
to n' = 500) is included in the exact cross section. Even though the approximation 
assumed only n' = 0, it does surprisingly well. For photon energies well above e^, the 
cross section tends toward the field-free relativistic (Klein-Nishina) cross section. The 
suppression of the relativistic magnetic cross section is due primarily to Klein-Nishina 
effects, i.e. electron recoil. The approximation for the differential scattering cross section 
for this case is given in Gonthier et al. (2000). 



5.4- Pair Production and Annihilation 

In a strong magnetic field, single photons as well as two or more photons may convert 
into electron-positron pairs. One-photon pair production cannot conserve both energy 
and momentum in field-free space, but a magnetic field can absorb the extra momentum 
of a photon with the energy required to created a pair. 



5.4-1 ■ One-photon pair creation and annihilation A photon with energy e, traveling 
at angle 9 to the magnetic field, can produce an electron with parallel momentum p 
and a positron with parallel momentum q only in the discrete Landau states that are 
kinematically allowed by the energy and momentum conservation equations, 

E„ + En' = e (90) 

p + q = e cos 6 (91) 

where En = {l+p'^ + 2nB'y^'^ and En' = (l+q"^ +2n' B'Y^'^ are the energies of the electron 
and positron. The threshold, e = 2/ sin 6', is the photon energy needed to produce a pair 
with momenta p = q = ecos9/2 in the ground state [n = n' = 0)||. In the case of high 
photon energies and low magnetic fields, e'^sm'^O/B' ^ 1, where the pair is produced 
far above threshold in high Landau states, the polarization-averaged pair production 
attenuation coefficient can be expressed in the asymptotic form (Klepikov 1954, Erber 
1966) 

mW = ^B'fAn9 T(x), (92) 

II This applies only to the || mode; for the _L mode, we must have n = 0, n' = 1 or n = 1. n' = 0, so 
the threshold is replaced by e sin 6=1 + \/l + 2B'. 
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Figure 6. Total polarization- averaged Compton scattering cross section for an electron 
at rest in the ground state, in units of the Thompson cross section, as a function of 
incident photon energy in units of cyclotron energy, for the case of incident photon 
angle sm9 = 0. Solid line: relativistic QED cross section (e.g. Bussard et al. 1986, 
Daugherty & Harding 1986) summed over final electron Landau states up to 500, 
dot- dashed line: magnetic Thompson cross section, Eqn. (|87|l . dotted line: Klein- 
Nishina (non-magnetic) cross section, dots: approximate QED cross section in Eqn 
(IHEJ- Adapted from Gonthier et al. (2000). 

f 0.377 exp (-7^) x«l 

T(x) ^ AMx-'/'Ai'ix-'^') = (93) 

[ 0.6x''/' X > 1 

where x = ^B' sin 0/2, Ai is the Airy function and A is the electron Compton wavelength. 
1 The probability of one-photon pair production thus rises exponentially with increasing 
photon energy and transverse field strength. A rule-of-thumb is that magnetic pair 
production will be important for this regime when the argument of the exponential in 
Eqn (|92p approaches unity, or when x ^ O-l- When the mag netic field S' ~ 0.1, this 
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condition will be satisfied below the threshold, so that in such high fields pair production 
will occur near threshold. The pair will then be produced in low Landau states and the 
attenuation coefficient will exhibit resonances at the threshold for each pair state (Toll 
1952, Daugherty & Harding 1983). Eqn will not be valid in this case, but may 
be corrected for near-threshold effects by making the substitution, x ~^ x/F-i where 
F = I + 0.42(a;sin6')"2-'^ in Eqn (Daugherty & Harding 1983), or by using the 
approximate expression (Baring 1988, see corrections in Baring 1991) 



Rl-{e, 9)^^ , [log (1/C)] 




1-1/2 



+ logC 



n-1/2 



4 



exp -i^-^7^1ogC , 2<ex«l/fi'. (94) 



where e_L = esin^, ^ = (e_L — 2)/(e_L + 2). Either of the above prescriptions will 
approximate the decrease in Ri^ near threshold (see Figure [7j). Semionova & Leahy 
(2001) have derived the one-photon pair production attenuation coefficient for all cases 
of photon polarization and electron and positron spin states, making use of the proper 
spin eigenstates defined by Sokolov & Ternov (1986). 

In the case where an electric field is present perpendicular to the magnetic field, 
the pair production attenuation coefficient can be obtained by a Lorentz transformation 
perpendicular to the magnetic field (Daugherty & Lerche 1975). When an electric field 
of strength E parallel to the magnetic field is present, Daugherty & Lerche (1976) 
determined that the pair attenuation coefficient is increased by an additional amount 
of order E'^/B'^, which is usually small in neutron star magnetospheres since E <^ B is 
always true for rotationally induced electric fields. 

The inverse process to one-photon pair creation is one-photon pair annihilation, 
in which an electron-positron pair annihilates into a single photon and is likewise 
only permitted in a strong field. The kinematic equations are the same as those for 
one-photon pair creation, Eqns ()90|) and (j^Tj) . with the electron and positron initially 
occupying Landau states n and n' with momenta p and q. The annihilation rate for pairs 
in the ground state, n = n' = has been calculated by Wunner (1979) and Daugherty 
& Bussard (1980) and the expression in the center-of-momentum frame (where p = —q) 
is particularly simple: 

RlTiE,) = ntn^ exp , (95) 

where A = 27rA is the electron Compton wavelength, Ep = [p^ + 1)^/^ and and n~ are 
the number densities of positrons and electrons. The annihilation rate in an arbitrary 
frame can be obtained by a Lorentz transformation along the magnetic field direction. 
Since the rate in Eqn increases exponentially with increasing field strength, one- 
photon annihilation overtakes the two-photon annihilation rate for pairs at rest at a field 
strength around 10^^ Gauss. One-photon annihilation from the ground state results in a 
line at 2m, broadened asymmetrically toward higher energies by the parallel momenta of 
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Figure 7. Onc-photon pair production rate as a function of plioton energy (in units 
of mc^), for two different magnetic field strengths (in units of the critical field). The 
solid lines are the approximation of Eqn H94|l to the exact rate and the dashed lines 
are the asymptotic formula of Eqn ()92|l . From Baring (1988). 

the pairs. Unlike in the case of two-photon annihilation (cf. ^5.4.2j) . Doppler broadening 
results only in a blueshift here, because the photon must take all of the kinetic energy of 
the pair in addition to the rest mass. The annihilation photons are emitted in a fan beam 
transverse to the field, which is broadened if the pairs have nonzero parallel momenta. 
Pairs annihilating from excited states produce additional lines above 1 MeV which at 
high energies blend together into a continuum. The one-photon annihilation rate of 
pairs from excited states can proceed at a rate orders of magnitude faster than from 
the ground state in fields below 10^^ G (Harding 1986, Wunner et al. 1986). The one- 
photon rate therefore becomes comparable to the two-photon rate (which does not rise 
rapidly in excited states) at lower field strengths for pairs in excited states. However, 
the effectiveness of synchrotron cooling may keep the densities of excited states low 
enough relative to the ground state to cancel the increase in the annihilation rate. The 
one-photon annihilation rate for different photon polarizations and particle spin states 
has been derived by Wunner et al. (1986) and Semionova & Leahy (2000). 

5.4-2. Two-photon pair creation and annihilation in a strong magnetic field Pair 
creation by two photons is significantly modified by strong magnetic fields from its 
field-free behavior. In field-free, the two-photon pair production cross section near 
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threshold in terms of the photon energy in the center-of-momentum frame, e^^ = 
[6162(1 - cos ei2)/2]V2^ is (Svensson 1982) 

'S'^X 121n(2e„) - l]/e«„, » 1 

where 61 and 62 refer to the energies of the photons and cos ^12 is the cosine of the angle 
between their propagation directions. The full relativistic cross section can be found in 
Jauch and Rohrlich (1980). In a magnetic field, the kinematic equations for this process 
impose conservation of energy and parallel momentum only: 

ei + 62 = K + En' (97) 

61 cos 9i + 62 cos 62 — p + q, (98) 

where 61 and 62 are their angles with respect to the field. The threshold depends on 
photon polarization direction with respect to the field, with the threshold for producing 
a pair in the ground state (n — n' — O), taking the form (Daugherty and Bussard 1980): 

(ei sin 9i + 62 sin 62)^ + 2e^e2 [1 - cos(^i - ^2)] > 4. (99) 

The second term is similar to the field-free threshold condition and the first term appears 
as a result of non conservation of perpendicular momentum. Thus it is possible for 
photons traveling parallel to each other (6*1 = 62 7^ 0°) to produce a pair, an event not 
permitted in field-free space. 

The two-photon pair production cross section in a strong magnetic field, like the 
one-photon pair production cross section, has resonances near threshold due to the 
discreteness of the pair states. The two-photon cross section in a strong magnetic field 
has been calculated by Kozlenkov and Mitrofanov (1987) for photon energies below one- 
photon pair production threshold and shows the same sawtooth behavior as the one- 
photon process. Near threshold, the magnetic field decreases the cross section below 
its free-space value, due to the decreased phase space available to the pair. Above its 
threshold, one-photon pair creation will dominate since it is a lower order process than 
two-photon pair creation. A comparison of the relative importance of the one-photon 
and two-photon processes (Burns and Harding 1984) shows that the one-photon process 
will generally dominate in magnetic fields above ~ 10^^ Gauss. 

The free-space cross section for two-photon pair annihilation in the non-relativistic 
limit is (727 ~ (3/8)ar//3r, where f3r is the relative velocity of the positron and electron. 
The annihilation rate for unpolarized positrons and electrons with densities n+ and n_ 
is then 

3 

-R27 = -crTcn_|_n_. (100) 
8 

The two-photon annihilation rate in a strong magnetic field has been calculated for 
pairs in the ground Landau state (Daugherty and Bussard 1980, Wunner 1979) and is 
unchanged from the field-free rate below 5' ~ 0.1. At 5' ~ 0.2, the rate decreases 
sharply due to the smaller phase-space of the virtual pair states, just as the one-photon 
annihilation rate is increasing exponentially. Two-photon annihilation of non-relativistic 
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Figure 8. Two-photon annihilation spectra for pairs at rest in a magnetic field in 
the ground Landau state, for different field strengths (in units of the critical field) and 
viewing angles to the field, d\. The kinematic cutoffs occur when | cos 02 1 = 1, where 
Bi is the angle of the unseen photon. From Baring & Harding (1992). 



pairs results in a line at 511 keV as in free space, but the relaxation of transverse 
momentum conservation in a magnetic field causes a broadening mostly on the red 
side of the line at viewing angles other than 90°. At a viewing angle of 90° to the 
field direction, the broadening is symmetric for i?'<0.1 but becomes asymmetrically 
broadened on toward the blue side for B' ~ 0.1, as shown in Figure |H1 The magnetic 
broadening can be approximated as 

572, sin^<V2^ 



(Ae) 



sin^./572, sin0>V25' 



(101) 



There is an increasing tendency in very high fields for one of the photons to be 
produced with almost all of the pair energy, so that two-photon annihilation behaves 
more like one-photon annihilation. The angular distribution of the photons from 
annihilation at rest also becomes more anisotropic with increasing field strength, with 
the peak of emission perpendicular to B , again similar to one-photon annihilation. 
Similar broadening occurs when the pairs have non-zero momenta parallel to the 
magnetic field (Kaminker et al. 1987, Baring & Harding 1992). 



5.5. Bound Pair Creation 

As was first noted by Shabad & Usov (1982), the dispersion curve (as shown in Figure 
in} of a photon propagating along a curved magnetic field, with e < 2/sin^ initially 
below the threshold for one-photon pair creation, will cross the lowest bound state of 
positronium prior to reaching the threshold for creating a free pair in the ground state 
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(n = n' ~ for || polarization, n = 0,n' = 1 or n = l,n' = for _L polarization). This 
situation was originally interpreted as a "capture" of the photon by the magnetic field 
for B' ~ 0.1, and a complete prevention of free pair creation until the photon could 
travel to a region of lower field strength. But Herold et al. (1985) pointed out that as 
the photon approaches the crossing point, it will adiabatically convert into positronium, 
after evolving through a mixed photon-positronium state (Shabad & Usov 1985, 1986), 
if the magnetic field is strong enough. The positronium state produced in this way from 
a single photon is stable to annihilation, since positronium cannot annihilate from the 
ground state in a magnetic field (Wunner et al. 1981). In this case, creation of a free 
pair is suppressed unless the positronium is unbound and several possible mechanisms 
have been discussed. Positronium may be ionized by an electric field parallel to the 
magnetic field, whose strength is sufficient when combined with the Coulomb field to 
form a potential barrier low enough for the electron and positron to tunnel through to 
become a free pair. The probability of ionization by a field of strength E\\ is (Usov & 
Melrose 1995) 



eEn 



-4(A£oo) 



3/2 



(102) 



_3(e^ll/mc2)A_ 

where A^qo is the positronium binding energy in the ground state. On the scale of a 
neutron star radius, R, the critical ionizing field strength, below which the ionization 
probability is negligible, is 

2 mc^(Aeoo)3/^ 
^11 ~3 iA 

and is a slowly varying function of B'. For B' ~ 0.2, e£'||°° ~ 10^"^ Volts cm"^. 
Alternatively, the bound pair could be photoionized by a UV radiation field of sufficient 
intensity (Hcrold et al. 1985). For example, in a blackbody radiation field of temperature 
Tg = T/10^ K the mean-free path for photoionization of positronium moving with 
Lorentz factor F in a magnetic field 5' ~ 0.1 is approximately (Bhatia et al. 1992) 

5x10^ cm (104) 

For T > 10^ K, the photoionization mean-free path of positronium will be short enough 
to free a bound pair in a fraction of a neutron star radius. 



5.6. Photon Splitting 

Photon splitting is a third-order process in which one photon divides into two or more 
photons in the presence of a strong external field. This process can be most easily 
understood in the context of the vacuum polarization effect on photon propagation in 
a strong magnetic field, which is due to the interaction of the photon virtual pairs 
with the field. To lowest order, photon splitting occurs when one member of a vacuum 
polarization pair radiates. While it is kinematically allowed in field-free space, it is 
forbidden by Furry's theorem (Furry 1937, Jauch & Rohrhch 1980), a QED charge 
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Figure 9. Dispersion curves for photons propagating in a strong magnetic field. 
The sloped straight line is the unperturbed dispersion curve (w/c)^ — ifj^ = K^. 
£n,n'{'nc,Pc) are the energies of the positronium states with quantum number ric and 
parallel momentum pc- The thin solid curves, labeled by nf, are the positronium 
dispersion curves for the lowest states having n ^ n' — and n + n' = 1 for the 
electron and positron Landau states, where ± are the parities. The hatched borders 
are the thresholds for production of a free pair in the two lowest states, which are the 
same as the positronium continua for each scries of states. The thick solid curves are 
the lowest dispersion curves of the mixed photon-positronium states with n = n' = 0. 
The thick dashed lines are the lowest dispersion curves for the n + n' = 1 mixed states. 
Prom Shabad & Usov (1985). 



conjugation symmetry stating that the S-matrix vanishes for any closed diagram with 
an odd number of vertices having only external photon lines. The symmetry is broken by 

a strong magnetic field. The photon splitting rate in a magnetic field was computed by 
Bialynicka-Birula & Bialynicki-Birula (1970), and Adler et al. (1970, 1971) in the weakly 
dispersive limit. In this limit, there are three photon polarization modes permitted 
by QED (i.e., they do not violate CP conservation): ||^-L|| and A 

simplified expression (for B'<1) for the polarization-dependent attenuation coefficient 
for photon splitting was derived by Baring (1991): 

RU, = iflr-xi, fl?-xx = 5 (i) ' 
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and 

where e is the incident photon energy propagating at angle 6 to the magnetic field of 
strength B'. Since there is no threshold for photon splitting the expression above is valid 
for arbitrarily small photon energies. For large values of > 1, the B' dependence in 
Eqn ()106|) saturates (i.e., the factor B'^ 1 for B' ^ 1). The full expression for the 
photon splitting attenuation coefficient is given by Adler (1971) and Stoneham (1979). 

Since photon splitting is a third order process, its rate is smaller than that of one- 
photon pair production above e = 2/ sin 6'. However, it can compete with one-photon 
pair production in neutron star magnetospheres, where photons are below pair creation 
threshold at their emission points (see As the photon travels in a region of curved 
field lines, its angle to the field increases and the photon may split before reaching pair 
creation threshold (Baring 1991, 1995, Harding et al. 1997). Since the photon splitting 
rate is very sensitive to field strength, its infiuence grows quickly with increasing field 
strength with the switch from pair creation to photon splitting being dependent on 
field geometry. For photons propagating in neutron star dipole magnetic field geometry, 
photon splitting becomes dominant as an attenuation mechanism for 5' ~ 1 (Baring & 
Harding 1998). 

Adler (1971) examined the kinematic selection rules for photon splitting in the 
weakly dispersive limit and found that only the -L— ^-HH mode conserved energy and 
momentum in the magnetized vacuum [another channel, _L— >-||_L, is also allowed 
kinematically, but the rate is suppressed by a factor of order n\\—l (see Eq.lUUj) compared 
to the -L-^llll channel, and thus is negligible in the weakly dispersive limit]. Usov 
(2002), including only the linear positronium contribution to the polarization tensor, 
has shown that Adler's kinematic selection rules hold in arbitrarily high magnetic fields. 
However, the moderate vacuum dispersion present in magnetar fields (see ^lOj) may 
critically depend on higher order (non-linear) contributions from the polarization tensor 
(e.g. Melrose and Parle 1983). If Adler's selection rules are in fact ubiquitous it has 
profound implications, since it means that the only polarization mode permitted to 
split produces only photons of the mode that cannot split. Such a situation will prevent 
pure photon splitting cascades (Baring 1995, Harding & Baring 1996), where the two 
product photons of an initial splitting split into two more photons, each of which split 
further, reprocessing the spectrum to lower and lower energies. Such a selection rule also 
leads to an extremely strong polarizing mechanism, since only one mode is attenuated 
and produces photon of the other mode. Bulik (1998) has derived the photon splitting 
rate including the effect of a plasma. He finds that the presence of a plasma changes 
the kinematic selection rules, allowing addition splitting modes to become important. 
However, these changes seem to operate only in a small region of phase space in neutron 
star atmospheres, between B' ~ 0.1 — 2, 6 < 30° and plasma density p > lgcm~^. 
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6. Physics of Neutron Star Interiors and Envelopes 
6.1. Overview 

Many observable properties of neutron stars, such as the mass-radius relation, the 
maximum rotation rate and the thermal evolution, depend on the properties of the 
star's interior (e.g., Baym & Pethick 1979; Shapiro & Teukolsky 1983; Glendenning 
2000; Heiselberg 2002; Lattimer & Parkash 2005). An important goal of neutron star 
astrophysics is to use various observations of neutron stars to probe the property of 
matter under extreme conditions, particularly at super-nuclear densities. 

A neutron star can be divided into several different regions: the atmo- 
sphere/envelope, the crust and the core. The atmosphere/envelope contains a negligible 
amount of mass, but plays an important role in determining the emergent photon spec- 
trum and flux (see The crust, constituting 1 — 2% of the stellar mass and extending 
~ 1 km below the surface, contains atomic nuclei embedded in an electron Fermi sea. 
Since a neutron star is formed as a collapsed, hot {ksT > 10 MeV) core of a massive 
star in a supernova explosion, the neutron star matter may be assumed to be fully cat- 
alyzed and in the absolute lowest energy state (e.g., Salpeter 1961; Baym, Pethick & 
Sutherland 1971). The dominant nuclei in the crust vary with density, and ranges from 
^^Fe at p ^ 8 X 10^ g cm^^ to neutron-rich nuclei with A ~ 200 and Z/A ~ 0.1 near the 
core-crust interface at p ~ 10^^ g cm~^ (see Douchin & Haensel 2001 for recent work 
on neutron star crust composition)^. This "standard" picture of neutron star crust, 
in which one particular kind of nucleus is present at a given density and these nuclei 
form a near-perfect crystal lattice, has been challenged by Jones (1999), who argued 
that thermodynamic fluctuations at the time of freezing would ensure the presence of 
several kinds of nuclei at a given pressure, yielding an impure solid. This would sub- 
stantially increase the resistivity of the crust, particularly at low temperatures, and thus 
greatly affect magnetic fleld evolution and heat conduction in the crust (Jones 2004). 
At densities above the 4 x 10^^ g cm~^, neutrons "drip" out of nuclei and permeate the 
crust (this region is called "inner crust"). At temperatures less than ~ 0.1 MeV, these 
neutrons are expected to form Cooper pairs and turn superfluid, largely decoupling dy- 
namically from the rest of the star. This has been invoked to explain pulsar glitches, 
quick increases of the observed rotation rate the star, which might be caused by a sud- 
den transfer of angular momentum from the more rapidly rotating neutron superfluid to 
the crust. Around the highest density region in the crust, the nuclei may assume rod- 
and plate-like shapes ("nuclear pasta"; Pethick & Ravenhall 1995), before an eventual 
transition to uniform nuclear matter of the liquid core. 

The liquid core (with densities ^ 2 x 10^'* g cm~^) consists of a mixture of 
neutrons, protons and electrons, which at progressively higher densities are joined by 
muons and other exotic particles such as strangeness-bearing hyperons and/or pion 

% Neutron stars that have undergone accretion are expected to have different crust compositions because 
of nuclear reactions and weak interactions during the accretion phase (e.g., Haensel & Zdunik 1990; 
Blaes et al. 1990; Schatz et al. 1999). 



Physics of Strongly Magnetized Neutron Stars 



43 



or kaon Bose condensates. Chemical equilibrium among these particles is established 
by weak interactions such as neutron beta decay [n — > p + e~ + u) and electron 
capture (e~ + p ^ n + u), and the nuclear symmetry energy plays an important role 
in determining the relative abundance of neutrons and protons. It is possible that both 
neutrons and protons in the core are in a superfluid/superconducting state, in which 
the neutron vorticity is concentrated in quantized vortex lines, and magnetic flux may 
be similarly concentrated in proton flux tubes. The transition temperatures for these 
superfiuid states are highly uncertain. Such a superfiuid would alter the specific heat 
and neutrino emissivities of the core, thereby affecting the thermal evolution of the 
neutron star (Yakovlov & Pethick 2004). 

At sufficiently high densities, free quarks may appear in the stellar core. The so- 
called quark stars may come in several forms depending on the details of the nuclear to 
quark matter phase transition (Weber 2005). A pure quark star (also called "strange" 
star) consists of up, down and strange quarks with electrons to fulfill charge neutrality; 
the star may be bare or may be enveloped in a thin nuclear crust (Alcock et al. 1986). 
A hybrid star has a core of quark matter and a mantle of nuclear matter (e.g., Alford 
et al. 2005). A mixed star has mixed phase of nuclear and deconfined quark matter 
over a range of density and radius (Glendenning 1992). It is likely that a condensate of 
quark Cooper pairs may appear at low temperatures characterized by a BCS gap usually 
referred to as color superconductivity (e.g., Alford et al 2001). Such a phase is predicted 
to have a unique property (e.g., the densest QCD phase, the so-called color-flavor locked 
phase, is a color superconductor but an electromagnetic insulator). 

An important probe of the neutron star interior is the mass-radius relation. Neutron 
star masses have been measured accurately from timing observations of binary radio 
pulsars, and less accurately from binary X-ray pulsars (Stairs 2004). Of particular 
interest is PSR J0751+1807 (with white dwarf companion in a 6 hour orbit), which 
has a mass of 2.1 ± O.2M0 (Nice et al. 2005). This measurement, together with that 
of Vela X-1 (1.86 ± O.I6M0), rules out a class of soft equations of state. The neutron 
star radius can be measured from observations of the star's surface emission, including 
possible spectral line identifications. This requires detailed modeling of neutron star 
atmospheres (®. Recent observations of quiescent neutron star binaries in globular 
clusters (with known distances) yield an effective radius Rod = -R(l + z) (where z is the 
gravitational redshift at the neutron star surface) in the range of 13 to 16 km (Rutledge 
et al. 2002; Heinke et al. 2003), consistent with many nuclear equations of state. 

The study of neutron star cooling can potentially provide useful information on 
the neutron star interior, such as superfluidity, direct vs. modified URCA processes 
(which depend on the symmetry energy of nuclear matter). Theoretical cooling curves 
can be compared with observations if the neutron star age can be estimated and the 
atmosphere emission spectrum can be properly interpreted (Yakovlev & Pethick 2004). 
Observations of pulsar glitches and possible precession also provide indirect constraints 
on the interior physics (e.g. Link 2003; Akgun et al. 2006). In the future, measurement 
of neutron-star moments of inertia (in the double pulsar system PSR J0737-3039; Lyne 
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et al. 2004) and detection of gravitational waves from coalescing neutron star binaries 
will lead to useful information of neutron star equation of state (e.g., Lai & Wiseman 
1996; Faber et al. 2002; Shibata et al. 2005). 

Because of its high density, the equation of state of the bulk interior of a neutron 
star is not strongly affected by the magnetic field (e.g., Broderick et al. 2000), unless 
one considers field strengths approaching the maximum "virial" value, 

^max ^ 10^^ ( ( ' G, (107) 



^1.4M0y VlOkm. 

which is set by equating the magnetic energy (47ri?V3)(5V87r) to the gravitational 
binding energy GM^/R. A strong magnetic field, however, can significantly influence 
the property of neutron star envelopes, which we discuss in the following. 

6.2. Free Electron Gas in Strong Magnetic Fields 

Before discussing various properties of neutron star envelopes, we summarize the basic 
thermodynamical properties of a free electron gas in strong magnetic fields at finite 
temperature T (see, e.g.. Landau & Lifshitz 1980 for the B = case). 

The number density Ue of electrons is related to the chemical potential /ig by 

"^=(2^„?/"Ioo^^^- « 

where gn is the spin degeneracy of the Landau level {go = 1 and gn = 2 for n > 1), and 
/ is the Fermi-Dirac distribution 



f 



1 + exp 



(109) 



with En given by Eq. (fTH|l . The electron pressure is given by 

Note that the electron pressure is isotropic ^. The grand thermodynamic potential is 
Q = —P^V, from which all other thermodynamic quantities can be obtained. Note that 
for nonrelativistic electrons (valid for Ep <^ nieC^ and <^ meC^), we use Eq. (jHJ for 
En, and the expressions for the density Ue and pressure Pe can be simplified to 

+ The transverse kinetic pressure Pei_ is given by an expression similar to Eq. Hll()|l . except that p^c^ is 
replaced by (p^c^) — n(3{meC^Y . Thus the kinetic pressure is anisotropic, with P^^ ~ ^ Pe±+MB, 
where M is the magnetization. When the electron gas is compressed perpendicular to B, work must be 
done against the Lorentz force density (V x ) x B involving the magnetization current. Thus there is a 
magnetic contribution to the perpendicular pressure, with the magnitude MB. The composite pressure 
tensor is therefore isotropic, in agreement with the thermodynamic result Pg = —D,/V (Blandford and 
Hernquist 1982). For a nonuniform magnetic field, the net force (per unit volume) on the stellar plasma 
is -VPe - V(BV87r) + (B • V)B/(47r). 
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where Xxe = {^Trh^ /meksTy/'^ is the thermal wavelength of the electron, and is the 
Fermi integral: 



X ' 



■ dx. 



(113) 



/o exp(x — y) + 1 

First consider degenerate electron gas at zero temperature. The Fermi energy 
(excluding the electron rest mass) Ep = fie{T = 0) — mgC^ = {meC^)eF is determined 
from 
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(114) 
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where Ag = h/{mf,c) is the electron Compton wavelength, (3 = B/Bq, and n^ax is set 
by the condition (1 + epY > (1 + 2nmax/5). The electron pressure is given by 

Xp^n) 



;i + 2n/5)V2 



where 



e(y) 



-yv/l + y2__ln(y + ^l + y2)^ 



(116) 



(117) 



which approaches |/^/3 for ?/ <^ 1. The critical "magnetic density" below which only 
the ground Landau level is populated (rimax = 0) is determined by (1 + ep^ = 1 + 2/3, 
which gives [see Eq. (E7jl ] 



Pb = 0.802Y^'^b^/'^ g cm-3 = 7.04 x 10^ Y^'^B%'^ g 



cm 



:ii8) 



where Ye = Z/A is the number of electrons per baryon. Similarly, the critical density 
below which only the n = 0,1 levels are occupied (rimax = 1) is Pbi = (2 + V^)pb = 
3.414 For p < ps, equation ()114|) simplifies to 



1/2 _ 

g cm . 
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p = 3.31 X 10^ Y^-^Bu [(1 + epY - 1 

For nonrelativistic electrons {ep -C 1), the Fermi temperature Tp 
{nieC^ /kB)tp is given by 

Tp = 2.70 B^i (FeP)' K (for p < pb), (120) 

where p is in units of 1 g cm^^. For p ^ pB, many Landau levels are filled by the 
electrons, Eqs. ()108|) and ()110|) reduce to the zero-field expressions. In this limit, the 
Fermi momentum p_p is given by 

Pf ^ 



Xp 

rUeC rUeC 

and the Fermi temperature is 



-{?>Ti^nef'^ = 1.009 X 10-2 (FeP) 



1/3 



{B = 0) (121) 
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l+a;|-l) ~3.0x lO^(reP)^^^ K, {B = 0) (122) 
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Figure 10. Temperature-density diagram illustrating the different regimes of 
magnetic field effects on the thermodynamic properties of a free electron gas. The 
solid lines are for B = 10^^ G, short-dashed lines for B = 10^^ G, and long-dashed 
lines for B = 10^** G. For each value of B, the vertical lines correspond to p = ps 
(the density below which only the ground Landau level is occupied by the degenerate 
electrons) and p = psi (the density below which only the n — 0,1 levels are occupied); 
the Fermi temperature is shown for p < ps and for ps < P ^ Pbi] the line marked 
by "T = Tb" [see Eq. H123|l ] corresponds to the temperature above which the Landau 
level effects are smeared out. The dotted line gives the Fermi temperature at i? = 0. 
The magnetic field is strongly quantizing when p ^ ps and T <, Tb , weakly quantizing 
when p ^ pb and T <, Tb , and non-quantizing when T ^ Tb or p 3> pb- 



where the second equahty apphes to nonrelativistic electrons [xp <C 1). Comparison 
between Eq. ()120|) and Eq. ()122|) shows that the magnetic field lifts the degeneracy of 
electrons even at relatively high density (see Fig. llOj) . 

Finite temperature tends to smear out Landau levels. Let the energy difference 
between the n = nmax level and the n = nmax + 1 level be AEb- We can define a 
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"magnetic temperature" 



Tb = -j^ = -j^ [^sjl + 2n^ax/5 + 2(3-^1 + 2n^,,/?j . (123) 

Clearly, Tp = Tb p = ps (see Fig. fTUI) . The effects due to Landau quantization are 
diminished when T > Tg. For p < pB, we have Tg = ( v^l + 2/3 — 1) {nieC^/kB), which 



reduces to Tb — huJce/kB for (3 = B / Bq <^ 1. For p ^ pb (or nmax ^ 1), equation 
becomes 




^ = 1.34 X 10« (1 + 4)"'/' K, (124) 



where m* = \Jm1 + {.Vf/cY = nie^Jl + Xp, with Zi;' given by Eq. ()121|) . 

There are three regimes characterizing the effects of Landau quantization on the 
thermodynamic properties of the electron gas (see Fig.^1 see also Yakovlev & Kaminker 
1994): 

(i) p ^ Pb and T <, Tb: In this regime, the electrons populate mostly the ground 
Landau level, and the magnetic field modifies essentially all the properties of the gas. 
The field is sometimes termed "strongly quantizing". For example, for degenerate, 
nonrelativistic electrons (p < pb and T <^Tp <^ nieC^/kB), the internal energy density 
and pressure are 

Ue = ^UeEp, (125) 
Pe = 2ue = "^n^Ep oc S-y. (126) 

These should be compared with the B = expression Pe = 2ue/S oc p^^^. Note that for 
nondegenerate electrons (T ^ Tp), the classical ideal gas equation of state, 

Pe = n.kBT, (127) 

still holds in this "strongly quantizing" regime, although other thermodynamic 
quantities are significantly modified by the magnetic field. 

(ii) p ^ Pb and T <, Tb'. In this regime, the electrons are degenerate (note that 
Tp > Tb when p > pB] see Fig. irUj) . and populate many Landau levels but the level 
spacing exceeds kBT. The magnetic field is termed "weakly quantizing". The bulk 
properties of the gas (e.g., pressure and chemical potential), which are determined by 
all the electrons in the Fermi only slightly affected by such magnetic fields. 
However, the quantities determined by thermal electrons near the Fermi surface show 
large oscillatory features as a function of density or magnetic field strength. These de 
Haas - van Alphen type oscillations arise as successive Landau levels are occupied with 
increasing density (or decreasing magnetic field). The oscillatory quantities are usually 
expressed as derivatives of the bulk quantities with respect to thermodynamic variables; 
examples include heat capacity, magnetization and magnetic susceptibility, adiabatic 
index {din P^/ din p), sound speed, and electron screening length of an electric charge in 
the plasma (e.g., Ashcroft & Mermin 1976; Blandford & Hernquist 1982; Lai & Shapiro 
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1991; Yakovlev & Kaminker 1994). With increasing T, the oscillations become weaker 
because of the thermal broadening of the Landau levels; when T ^ Tb, the oscillations 
are entirely smeared out, and the field- free results are recovered. 

(iii) T ^ Tb OT p ^ pb'- In this regime, many Landau levels are populated and 
the thermal widths of the Landau levels (~ /c^T) are higher than the level spacing. 
The magnetic field is termed "non-quantizing" and does not affect the thermodynamic 
properties of the gas. 

6.3. Magnetized Crusts, Effects on Thermal Structure and Cooling of Neutron Stars 

As discussed above, the effect of Landau quantization on the density-pressure relation 
is important only for p ^ pb [see Eq. flllSj) ]. Deeper in the neutron star envelope 
(and the interior), we expect the magnetic field effect on the bulk equation of state 
to become negligible as more Landau levels are filled. In general, we can use the 
condition pB ^ p, or B12 ^ 27 (FePe)^^^, to estimate the critical value of B above which 
Landau quantization will affect physics at density p. For example, at B ^ 10^^ G, the 
neutronization transition from ^^Fe to ^^Ni (at p = 8.1 x 10^ g cm~^ for B = 0; Baym, 
Pethick & Sutherland 1971) in the crust can be significantly affected by the magnetic 
field (Lai & Shapiro 1991). 

The ions in the neutron star envelope form a one-component plasma and are 
characterized by the Coulomb coupling parameter 



where = (3/47™j)^/'^ is the Wigner-Seitz cell radius, Ui = p/nii = p/{Amp) is the ion 
number density, pg = p/(10^g cm~^) and Tg = T/(10^K). For F <^ 1, the ions form a 
classical Boltzmann gas whose thermodynamic property is unaffected by the magnetic 
field. For F ^ 1, the ions constitute a strongly coupled Coulomb liquid. The liquid 
freezes into a Coulomb crystal at F = F^ — 175, corresponding to the classical melting 
temperature Tm (e.g., Nagara et al. 1987; Potekhin & Chabrier 2000). The quantum 
effects of ion motions (zero-point vibrations) tends to increase F^ (Chabrier et al. 1992; 
Chabrier 1993) or even suppress freezing (e.g., Ceperley & Alder 1980; Jones & Ceperley 
1996). At zero- field, the ion zero-point vibrations have characteristic frequency of order 
the ion plasma frequency Qp, with 



For T -C Toebye ~ hupi/kB, the ion vibrations are quantized. The effects of magnetic 
field on strongly coupled Coulomb liquids and crystals have not been systematically 
studied (but see Usov et al. 1980). The cyclotron frequency of the ion is given by 
huci = h{ZeB / Ampc) = 6.3 (Z/A) B12 eV. The ion vibration frequency in a magnetic 
field may be estimated as (Up^ + cj^J^/^. Using Lindeman's rule, we obtain a modified 




(128) 




(129) 
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melting criterion: 



For ujci ^ ^pi, or B12 ^100 pg 




~ 175. 



, the magnetic field does not affect the melting criterion 



(130) 



and other properties of ion vibrations. 

A strong magnetic field can significantly affect the transport properties and thermal 
structure of a neutron star crust. Even in the regime where the magnetic quantization 
effects are small (p ^ ps), the magnetic field can still greatly modify the transport 
coefficients (e.g., electric conductivity and heat conductivity). This occurs when the 
effective gyro- frequency of the electron, u*^ = eB/{mlc), where m* = ^m"^ + {pp/cy, 
is much larger than the electron collision frequency, i.e.. 



where tq is the effective electron relaxation time. In a degenerate Coulomb plasma with 
a nonquantizing mag netic field, r ~ Sirh^ /{AZmle^ In A) ~ 6 x 10-i^me/(m*Zln A) s. 
When u!*^tq ^ 1, the electron heat conductivity perpendicular to the magnetic field, 
K_L, is suppressed by a factor {u*^tq)~'^. In this classical regime, the heat conductivity 
along the field, k\\, is the same as the B = value. In a quantizing magnetic field, 
the conductivity exhibits oscillatory behavior of the de Haas - van Alphen type. On 
average, the longitudinal conductivity is enhanced relative to the B = value. The 
most detailed calculations of the electron transport coefficients of magnetized neutron 
star envelopes are due to Potekhin (1999), where earlier references can be found (see 
also Hernquist 1984, Yakovlev & Kaminker 1994). 

The thermal structure of a magnetized neutron star envelope has been studied by 
many authors. Hernquist (1985) and Van Riper (1988) considered the cases where B 
is normal or tangential to the stellar surfaces, while Schaaf (1990), Heyl & Hernquist 
(1998,2001), Potekhin & Yakovlev (2001) (for a Fe crust) and Potekhin et al. (2003) 
(for an accreted magnetized envelope) analysed arbitrary magnetic orientations and 
used increasingly more accurate transport coefficients. In general, a normal magnetic 
field reduces the thermal insulation as a result of the (on average) increased k\\ due to 
Landau quantization of electron motion, while a tangential magnetic field (parallel to 
the stellar surface) increases the thermal insulation of the envelope because the Larmor 
rotation of the electron significantly reduces the transverse thermal conductivity k,±. A 
consequence of the anisotropic heat transport is that for a given internal temperature 
of the neutron star, the surface temperature is nonuniform, with the magnetic poles 
hotter and the magnetic equator cooler (see Fig. [TT|) . Additional heating and surface 
temperature inhomogeneity comes from the bombardment of high-energy particles from 
the magnetosphere on the polar cap. Models for particle acceleration in neutron star 
magnetospheres predict that some fraction of electron-positron pairs that are produced 
in pair cascades (see ^ decelerate, turn around and accelerate downward to the 
neutron-star surface. In space-charge limited flow models (e.g. Arons & Scharlemann 
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Figure 11. Temperature profiles througli tlie iron envelope of a neutron star at (from 
left to right) B = lO", 10^^ lO" and 10^^ G. The internal temperature is fixed to 
10^ K (solid lines) or 10^ K (dot-dashed lines). The lines of each group correspond to 
cosO — 1 (the lowest line), 0.7, 0.4, 0.1, and (the highest line). From Potekhin & 
Yakovlev (2001). 



1979, Harding & Muslimov 2001), the trapped positrons needed to screen the E\\ are 
only a small fraction of the number of primary electrons. Such models predict moderate 
heating, so that the heated polar caps have X-ray luminosities that are small compared 
to cooling luminosities for typical ages less than about 10^ yr. In polar vacuum gap 
(e.g. Ruderman & Sutherland 1975) and outer gap (Cheng et al. 1986) models, roughly 
half of the produced particles return to heat the surface, causing a much higher degree 
of polar cap heating. 

A superstrong magnetic field {B > 10^"^ G) also affects the cooling curve of a 
neutron star (see Potekhin & Yakovlev 2001; Potekhin et al. 2003). This is because for 
a given core temperature (whose time evolution largely depends on neutrino emission 
from the core, and thus is unaffected by the field strength unless B ^ 10^^ G), an 
enhanced radiation flux emerges from the magnetic polar regions, and such enhancement 
more than compensates for the reduced flux from the equatorial region where the heat 
conductivity is suppressed. Figure IT^ shows some examples of theoretical cooling curves 
for neutron stars with mass 1.5M0 and 1.3Mq, and with different surface field strengths 
(assuming dipole field structure). Note that in these models, the 1.5M0 star has a 
sufficiently high central density for the direct URCA process to operate, while only the 
slower, modified URCA process operates in the l.SM© star (see Yakovlev & Pethick 
2004); thus the 1.5M0 star cools much faster than the 1.3Mq star. We also see that 
for t < 10^ years, a superstrong magnetic field {B > 10^^ G) makes the neutron star 
brighter, while at later times the star cools down more quickly compared to the B = 
case. 




Figure 12. Effective surface temperature (as seen by distant observer) versus neutron 
star age for assumed mass 1.3Mq and 1.5Mq. The dots with error bars show estimates 
of neutron star ages and effective temperatures from various observations; the dots 
with errors indicate observational upper hmits. The different curves show the cooling 
of neutron stars with iron envelope for different magnetic field strengths (log B in 
Gauss). From Chabrier, Saumon & Potckhin (2006). 



6.4- Outermost Envelopes 

The outermost layer of a neutron star is of great importance since it mediates the 
emergent radiation from the stellar surface to the observer. The chemical composition 
of this layer is unknown. The surface of a fully catalyzed neutron star is expected to 
consist of ^^Fe formed at the star's birth, although C, O may also be present. This 
may be the case for young radio pulsars that have not accreted any gas. Once the 
neutron star accretes material, or has gone through a phase of accretion, either during 
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supernova fallback, or from the interstellar medium or from a binary companion, the 
surface composition can be quite different. A H layer is expected even if a small amount 
of fallback/accretion occurs after neutron star formation. A He layer results if H is 
completely burnt out. The atmosphere composition may also be affected by the (slow) 
nuclear burning on the NS surface layer (Chang, Arras & Bildsten 2004) and by the 
bombardment of high-energy particles on the surface (Thompson & Beloborodov 2005). 

Despite the uncertainties in the surface composition mentioned in the last 
paragraph, a great simplification arises due to the efficient gravitational separation 
of light and heavy elements. As an example, consider a trace amount of ions (with 
charge number Zi and mass number Ai) embedded in the ionized hydrogen plasma 
(with the electron number density rig) of a neutron star atmosphere. The downward 
drift speed of the ion relative to the H plasma is Wdrift — gtf, where g is the gravitational 
acceleration. The friction timescale tf can be estimated from AiTUp/tf ~ UerUpapiVp, 
where the ion-proton Coulomb collision cross-section is api ~ {ZiC^ /ksTY In A (here In A 
is the Coulomb logarithm), and the thermal velocity of the proton is Vp ~ {ksT /mpY^"^ . 



Using typical parameters ((7 ~ 2 x 10^^ cm s~^, T ~ 10^ K, ~ 10^^ cm~^), we find 
that the drift speed is of order 10 {A-J Z'fhiA) cm s^^. The timescale for gravitational 
settling over the pressure scale height H ~ ksT / {nipg) is of order a second. Thus the 
lightest elements, H and He, if present, are the most important chemical species in the 
outermost layer of the star. 

The hydrogen envelope of a neutron star may assume different forms, depending on 
the temperature and magnetic field. For "modest" field strengths {B <, a. few x 10^^ G) 
and temperatures (T ^ 10^ K), conditions satisfied by most observable neutron stars, the 
atmosphere is largely nondegenerate, and consists mainly of ionized hydrogen, H atoms 
and small H^v molecules, and the condensed phase can be neglected in the photosphere. 
As density increases, the matter gradually transform into a degenerate Coulomb plasma 
(see above). The atmosphere (with density p ~ 10~^ — 10^ g cm~^) constitutes a non- 
ideal, partially ionized. Coulomb plasma. An important issue concerns the ionization 
equilibrium of atoms in such a plasma. This is nontrivial because in a strong magnetic 
field, an atom cannot move free across the field, and there is strong coupling between 
the internal atomic structure and the center-of-mass motion (see ^3.H) . Lai & Salpeter 
(1995) gave an approximate analytic solution for a limited magnetic- field-temperature- 
denisty regime (see also Khersonskii 1987; Pavlov & Meszaros 1993). To date the most 
complete treatment of the equation of state (including Saha equilibrium) of partially 
ionized hydrogen plasma is that of Potekhin et al. (1999) (for 10^^ B <^ 10^"^ G) and 
Potekhin & Chabrier (2004) (for 10^'^ G^ B <^ 10^^ G). They used the numerical energy 
levels and fitting formulae of a moving H atom as obtained by Potekhin (1994,1998), 
and their free energy model takes into account the Coulomb plasma nonideality, the 
interactions of bound species with one another and with the electrons and protons. 



This gives 



Z^lnA 



Ai ni. 
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Figure 13. Domains of partial ionization in the p — T plane for B — 10^^, 10^'^ and 
10^** G. The contours delimit the domains where the atomic fraction xh < 0.1% (to 
the right of the dotted lines), 0.1%< xh < 1% (between the dashed and dotted lines), 
1%< xh < 10% (between the solid and dashed lines) or xh > 10% (to the left of the 
solid lines). From Potekhin et al. (2006). 



Ionization equilibrium is given by minimization of the free energy with respect to particle 
numbers under the stoichiometric constraints (keeping the volume and the total number 
of free and bound protons constant). The derivatives of the free energy with respect 
to p and T and their combinations provide the other thermodynamic functions. Figure 
El shows the domains of partial ionization in the p — T plane for different values of B. 
The higher the greater T at which the bound species are important. The calculated 
atomic fractions are needed to obtain the radiative opacities in the atmosphere. Note 
that the Potekhin et al. model becomes less reliable at relatively low T where larger 
molecules can be important (Lai 2001). Generally, this occurs within the p — T domain 
where xr ^ 0.1 (i.e., to the left of the solid lines in Fig. HSl). 

For sufficiently strong magnetic fields and/or low temperatures, the hydrogen layer 
may be in a condensed, metallic form, since the binding energy of the condensed 
hydrogen increases as a power-law function of -B, while the binding energies of atoms 
and small molecules increase only logarithmically. Lai & Salpeter (1997) studied the 
phase diagram of H under different conditions and showed that in strong magnetic fields, 
there exists a critical temperature Tcrit below which a phase transition from gaseous to 
condensed state occurs, with kTcrit about 10% of the cohesive energy of the condensed 
hydrogen. Thus, Tc^t ~ 8 x 10^ 5 x 10^ 10^ K for 5 = lO^^^ 10l^ 5 x 10^^ G (Lai 
2001). An analogous "plasma phase transition" was also obtained in an alternative 
thermodynamic model for magnetized hydrogen plasma (Potekhin et al. 1999). While 
this model is more restricted than Lai & Salpeter (1997) in that it does not include Hjv 
molecules, it treats more rigorously atomic motion across the strong B field and Coulomb 
plasma nonideality. In the Potekhin et al. model, the density of phase separation is 
roughly the same as in Lai & Salpeter (1997), i.e., ps^ ~ 1^^ B12 g cm"^ [see Eq. (jSH)], 
but the critical temperature is several times higher. Thus, although there is a factor of 
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a few uncertainty in T^rit, it is almost certain that for T < Tcrit/'2, the H surface of a 
neutron star is in the form of the condensed metalhc state, with neghgible vapor above 
it. The radiative property of such condensed phase is discussed in §H1 

For heavy elements such as Fe, detailed study of the phase diagram is not available 
since our current knowledge of the various states/phases in strong magnetic fields is not 
as complete as for H (see ^3.4|) . Calculations so far have shown that at B = 10^^ — 10^^ G, 
a linear chain is unbound relative to individual atoms for Z ^ 6, and thus chain-chain 
interactions play an important role in determining whether 3D zero-pressure condensed 
matter is bound or not. For larger field strengths, B ~ 10^^(Z/26)^ G, a linear chain 
becomes bound (Medin & Lai 2006). If the condensed Fe is unbound with respect to 
the Fe atom, then the outermost Fe layer of the neutron star is characterized by gradual 
transformation from nondegenerate gas at low densities, which includes Fe atoms and 
ions, to degenerate plasma as the pressure (or column density) increases. On the other 
hand, even a weak cohesion of the Fe condensate can give rise to a phase transition 
at sufficiently low temperatures. Crude numerical results by Jones (1986) together 
with approximate scaling relations suggest an upper limit of the cohesive energy (for 
Z > 10) Qs ^ Z'^/^b12 eV. Thus for Fe, the critical temperature for phase transition 



The iron surface layers of magnetic neutron stars can also studied using Thomas- 
Fermi type models (e.g., Fushiki, Gudmundsson & Pethick 1989; Thorolfsson et 
al. 1998). While these models are too crude to determine the cohesive energy of 
the condensed matter, they provide a useful approximation to the gross properties of 
the neutron star surface layer. At zero temperature, the pressure is zero at a finite 
density which increases with increasing B. This feature is qualitatively the same as 
in the uniform electron gas model. Neglecting the exchange- correlation energy and 
the nonuniformity correction, we can write the pressure of a zero-temperature uniform 
electron g as as 



where the first term is given by Eq. pi6|l [or by Eq. p26j) in the strong field, degenerate 
limit], and the second term results from the Couloumb interactions among the electrons 
and ions. Setting P = gives the condensation density ps,o as in Eq. (jHT|) . At finite 
temperatures, the pressure does not go to zero until p — > 0, i.e., an atmosphere is 
present. Note that the finite-temperature Thomas- Fermi model only gives a qualitative 
description of the dense atmsophere; important features such as atomic states and 
ionizations are not captured in such a model. 

7. Neutron Star Magnetic Field Evolution 

Magnetic field is perhaps the single most important quantity that determines the various 
observational manifestations of neutron stars. Thus it is natural that a large amount 
of work has been devoted to the study of neutron star magnetic field evolution. Recent 
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reviews include Bhattacharya & Srinivasan (1995), Ruderman (2004), Reisenegger et 
al. (2005). 

7.1. Observations 

As discussed in ^ with a few exceptions, our current knowledge of neutron star 
magnetic fields are largely based on indirect inferences: 

(i) For radio pulsars, of which about 1600 are known today (see Manchester 2004), 
the period P and period derivative P can be measured by timing the arrivals of radio 
pulses. If we assume that the pulsar spindown is due to magnetic dipole radiation, we 
obtain an estimate of the neutron star surface field B [see Eq. (^]. For most radio 
pulsars, the inferred B lies in the range of 10^^ — 10^^ G. For a smaller population of 
older, millisecond pulsars, B ~ 10^^^ G - Such field reduction is thought to be related 
to the recycling process that turns a dead pulsar into an active, millisecond pulsar. 

(ii) For Anomalous X-ray pulsars and Soft Gamma Repeaters (Woods & Thompson 
2005), superstrong magnetic fields {B ~ lO^'' — 10^^ G) are again inferred from 
the measured P (in the range 5 — 12 s) and P (based on X-ray timing study) and 
the assumption that the spin-down is due to the usual magnetic dipole radiation or 
Alfven wave emission. In addition, strong theoretical arguments have been put forward 
(Thompson & Duncan 1995,1996,2001) which suggest that a superstrong magnetic field 
is needed to explain various observed properties of AXPs and SGRs (e.g., energetics of 
SGR bursts/flares, including the spectacular giant flares from three SGRs: the March 
5 1979 flare of SGR 0525-66 with energy E ^ Q x 10^^ erg, the August 27 1998 flare 
from SGR 1900+14 with E^2x 10^^ erg, and the December 27 2004 flare from SGR 
1806-20 with ~ 4 x 10^^ erg; quiescent luminosity ~ IQ'^^"^^ erg/s, much larger 
than the spindown luminosity IflCl) (see ^10|). Tentative detections of spectral features 
during SGR/AXP bursts have been reported in several systems (e.g., Gavriil et al. 2003; 
Ibrahim et al. 2003; Rea et al. 2003), which, when interpreted as proton cyclotron lines, 
imply B ~ 10^^ G. 

(iii) For a number of thermally emitting neutron stars ("dim isolated neutron 
stars"), spectral features with energy 0.2 — 1 keV have been detected (see §8.1|) . While 
the identification of these features is unclear, if one assumes that they are due to 
proton cyclotron resonance or atomic transitions in light elements, one would infer 
B ~ 10^3 - 10^4 G. 

(iv) For about a dozen accreting X-ray pulsars in binary systems, electron cyclotron 
features have been detected, implying B ~ 10^^ — 10^^ G (see §1). For many other X- 
ray pulsars with no detectable cyclotron features, we can use the measured spin period 
together with theoretical ideas of spin equilibrium (i.e., spinup due to accretion of matter 
is balanced by magnetic braking - loss of stellar angular momentum via magnetic fields) 
to estimate the surface magnetic field of the neutron star. For the ~ 100 neutron stars in 
high-mass X-ray binaries, this typically gives B ~ 10^^ G (e.g., Bildsten et al. 1997). In 
the last few years, half a dozen or so accreting millisecond pulsars have been discovered 
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(e.g., Chakrabarty 2005), with the estimated field strength B ~ 10^ G — these systems 
are thought to be neutron stars undergoing the rccychng process. 

Several lines of observations/arguments point toward the possibility of an evolving 
magnetic fields in neutron stars. 

(i) From the discussion above, it appears that young neutron stars have strong 
magnetic fields ~ 10^^ — 10^^ G (most radio pulsars, magnetars, high-mass X-ray 
binaries), whereas old neutron stars have weak fields ^ 10^ G (millisecond pulsars, 
low-mass X-ray binaries). If these two groups have an evolutionary connection, then 
the field must be significantly reduced on timescale of 10^ yr. Most hkely, this reduction 
is associated with the accretion process that recycled the neutron star. 

(ii) Magnetic field evolution in isolated pulsars has long been studied using the 
statistical distribution of pulsars on the P — P diagram, as well as pulsars' spatial 
and velocity distributions. Such study probes field evolution on the timescale of the 
ages of radio pulsars (^ 10^ yr). Unfortunately, the conclusions of such population 
study have often been conflicting. For example, Narayan & Ostriker (1990) and 
Gonthier et al. (2004) suggested that field decay occurs on timescale ^ 5 Myr. Similar 
studies by Bhattacharya et al. (1992), Lorimer et al. (1997), Faucher-Giguere & Kaspi 
(2006) found no evidence for field decay during the radio pulsar lifetime, implying a 
decay time constant ^ 100 Myr. Such divergent results illustrate the difficulties of 
controlling various systematic uncertainties (e.g. selection effects, luminosity evolution 
law, dependence of beaming fraction on period) in pulsar population studies. 

(iii) Magnctar emission is most likely powered by magnetic energy, implying the 
dissipation of superstrong magnetic fields (see Thompson & Duncan 1995,1996; Woods 
& Thompson 2005). 

(iv) One explanation for the "anomalous" braking indices n < 3 in young pulsars 
is that their magnetic dipole field increases with time (but see Melatos 1997). It has 
also been pointed out that following each pulsar glitch, the apparent magnetic field 
cx: ^/PP) increases, which might imply that some pulsars (Crab and Vela), may evolve 
into magnetars (Lync 2004; Lin & Zhang 2004). Since many AXPs and SGRs are 
observed to be associated with supernova remnants or clusters of massive stars, most 
magnetars are expected to be younger than radio pulsars, rather than the other way 
around. 

7. 2. Origin of Neutron Star Magnetic Field 

The magnetic fields of neutron stars were most likely already present at birth. The 
traditional fossil field hypothesis suggests that the magnetic field is inherited from the 
progenitor, with magnetic fiux conserved and field amplified {B oc R"^) during core 
collapse. In the case of magnetic white dwarfs (with measured fields in the range 
~ 3 X 10^ — 10^ G), there is strong evidence that the fields are the remnants from 
a main-sequence phase (Ap/Bp stars, with B ~ 200 G -25 kG) (see, e.g., Ferrario 
& Wichramasinghe 2005). Neutron stars descend from main-sequence stars with mass 
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> 8Mq (i.e., O and early B stars). Only recently have the large-scale magnetic fields 
(with B 1 kG) of O stars been detected (in two stars so far; see Donati et al. 2006). 
It is interesting that the magnetic flux of such O stars $ ~ IO^ttRq G (for R ~ lOi?©; 
of course, not all of this fiux threads the inner 1.4Mq core) is of the same order of 
magnitude as the flux of a 10^^ G neutron star {R ~ 10~^Rq) as well as the fluxes 
of the most strongly magnetic Ap/Bp stars and white dwarfs. It is also of interest to 
note that magnetic white dwarfs {B ^ 1 MG) tend to be more massive (mean mass 
~ O.93M0) than their non- magnetic counterparts (mean mass ~ O.6M0). Since white 
dwarfs with M ^ O.7M0 and neutron stars form exclusively from the material that 
belongs to the convective core of a main-sequence star (Reisenegger 2001), this suggests 
that the magnetic field may be generated in the convective core of the main-sequence 
progenitor. It has also been suggested that magnetized neutron stars could be produced, 
in principle, by accretion- induced collapse of magnetic white dwarfs (Usov 1992). 

Alternatively, it has been argued that magnetic field may be generated by a 
convective dynamo in the first ~ 10 seconds of a proto-neutron star (Thompson & 
Duncan 1993). In principle, the maximum field achievable is either B ~ (47rp)^/^^;con ~ 
4 X 10^^ G (for convective eddy speed Vcon ~ 10^ km s~^) or S ~ {Anpy^'^RAQ ~ 
2 X 10^^ G [for differential rotation AQ ~ 27r/(l ms)]. How much dipole field can be 
generated is more uncertain. It could be that a large-scale field of ~ 10^^ G is generated 
if the initial spin period of the neutron star is comparable to the convective turnover 
time (~ Hp/ Vcon ~ 1 ms for pressure scale height of 1 km), whereas in general only 
a small-scale (~ Hp) field is produced, resulting in mean field of order 10^^ — 10^^ G 
typical of radio pulsars. We note that studies of presupernova evolution of massive stars, 
including magnetic fields, typically find that a newly formed neutron star has rotation 
rate appreciably slower than the breakup rate (Heger et al. 2005). 

In either the fossil field or dynamo scenario, the magnetic field of a proto-neutron 
star is likely "messy" and unstable. Before the crust forms (around 100 s afetr collapse), 
the field will evolve on the Alfven crossing time tA ~ {A'kpY'^R/B (~ 0.1 s for 
B ~ 10^^ G) into a stable configuration in which a poloidal field coexists with a toroidal 
field of the same order of magnitude (see Braithwaite & Spruit 2006). 

Gradual field generation due to a thermomagnetic effect in the crust has been 
discussed (Urpin & Yakovlev 1980; Blandford et al. 1983; Wiebicke & Geppert 1996), 
but it seems unlikely to be capable of producing large-scale fields stronger than 10^^ G. 

7.3. Physics of Field Evolution in Isolated Neutron Stars 

The magnetic field in a neutron star evolves through a series of quasi-equilibrium states, 
punctuated by the release of elastic stress in the crust and hydro dynamic motions in 
its liquid core. The physics of the quasi-equilibrium field evolution was discussed by 
Goldreich & Reisenegger (1992) (see also Reisenegger et al. 2005). The bulk region 
of a neutron star comprises a liquid core of mostly neutrons with a small fraction 
(Fg ~ a few%) of protons and electrons. The medium is stably stratified due to the 
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Yg gradient, thus the magnetic field cannot force the bulk neutron fiuid to move (e.g. 
due to buoyancy) unless -B^/(87r) ^ Y^P (or B ^ 10^^ G) or the fiuid can change its 
composition as it rises (which takes place on a timescale longer than the neutrino cooling 
time for B ^ 10^^ G). The magnetic field evolves according to 

f = -Vx(^VxB)+Vx(--i-.B)+Vx(v.xB). (134) 

where the three terms on the right-hand side represent three different effects: 

(i) The first term represents Ohmic diffusion of the magnetic field (cr is the zero-field 
conductivity defined by j = aE) , with timescale 

^Ohmic ~ n — (135j 

for fields which vary on lengthscale L. In the core, the conductivity is large since all 
the particles are degenerate, leading to tohmic ~ 2 x 10^^{L]^in/Tsy{p/ Pnuc)^ yr (Baym, 
Pethick & Pines 1969), where pnuc = 2.8 x 10^'' g cm'^, Lkm = L/{1 km). Of course, 
currents which are confined to the crust would have a much shorter decay time (e.g.. 
Sang & Chanmugam 1987; Gumming et al. 2004). 

(ii) The second term describes advection of the field by Hall drift: The magnetic 
field is carried by the electron fiuid, which drifts with respect to the ions with velocity 
Ve = — j/(nee). This term is non-dissipative, but can change the field structure on 
timescale 

t„.,~^^~5xlO»|H( M (136) 

CB Bi2 VPnuc/ 

Goldreich & Reisenegger (1992) suggested that the nonlinear Hall term may give rise 
to a turbulent cascade to small scale, thus enhancing the Ohmic dissipation rate of the 
field. This "Hall cascade" has been confirmed by recent simulations of electron MHD 
turbulence (also known as whistler turbulence; see Gho & Lazarian 2004). Several 
studies have also demonstrated magnetic energy transfer between different scales (e.g., 
Urpin & Shalybkov 1999; Hollerbach & Riidiger 2002). A specific example of the Hall- 
drift enhanced Ohmic dissipation can be seen by considering a purely toroidal field (see 
Reisenegger et al. 2005): In this case, the evolution equation for the field reduces to 
the Burgers equation (see also Vainshtein et al. 2000), where the field tends to develop 
a cusp (and the associated current sheet) and fast Ohmic decay. Recently, Rheinhardt 
& Geppert (2002) and Rheinhardt et al. (2004) discussed a "Hall drift instabihty" that 
can lead to nonlocal transfer of magnetic energy to small scales, leading to enhanced 
crustal field dissipation. 

(iii) The third term describes ambipolar diffusion, which involves a drift of the 
combined magnetic field and the bulk electron-proton fiuid relative to the neutrons. 
The drift velocity is determined by force balance m.pVa/'Tpn = ~ V(Ayu), where 
Tpn is the proton-neutron collision time, = j x B/(cnp) is the magnetic force per 
proton-electron pair, and V(A/i) (with A/x = /Xp + /ig ~ A'-n) is the net pressure force due 
to imbalance of /3-equihbrium (Goldreich & Reisenegger 1992). There are two modes 
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of ambipolar diffusion: Tlie solenoidal mode (V ■ = 0) is noncompressive, and does 
not perturb /3-equilibrium, thus = Tp^Tp^fs/mp, and the associated diffusion time 

4mb ~ L/Va is 

tU-^^^^^-3xlO^S^yr, (137) 

where we have used Up ~ 0.05p/mp and rp„ ~ 2x 10^^''T8"^(p/pnuc)^^^ s. The irrotational 
mode (V X Vq = 0) is compressive, and is impeded by chemical potential gradients 
V(A/i), and is possible only if weak interaction re-establishes /9-equilibrium during the 
drift. This gives a timescale t^mb ~ icooi/ Bf^ ~ 5 x 10^^T^^Bi2 yr (where tcooi is the 
neutrino cooling time and modified URCA rate has been used). Like Ohmic decay, 
ambipolar diffusion is dissipative, which leads to heating of the core and deep crust of a 
magnetar; this may be the power source for the persistent X-ray emission of magnetars 
(Thompson & Duncan 1996; Heyl & Kulkarni 1998; Arras et al. 2004; see 

In addition the the "steady" field evolution discussed above, crust fracture can 
also lead to sudden change of the magnetic field (Thompson & Duncan 1996). The 
crust has a finite shear modulus fi, and when the yield strain 6'max is exceeded the 
lattice will fracture. The characteristic yield field strength is of order ~ (47r^max/^)^^^ = 
2 X 10^^(^max/10~^)^/^ G. For example. Hall drift of the magnetic field causes the stresses 
in the crust to build up, and irregularities in the field can be damped quickly by crustal 
yielding - this may be responsible for magnetar bursts (Thompson & Duncan 1995; see 

m- 



7.4- Accretion- Induced Field Reduction 

The weak (10^ — 10^ G) magnetic fields found in millisecond pulsars and low-mass X-ray 
binaries have long led to the suggestion that the dipole field is reduced during accretion, 
perhaps through diamagnetic screening of the magnetic flux by the accreted plasma 
(e.g., Bisnovati-Kogan & Komberg 1975; Shibazaki et al. 1989; Romani 1990). Exactly 
how this is achieved (if at all) is not clear, and a number of idealized problems have been 
studied. To give some recent examples: Gumming et al. (2001) considered a ID model in 
which matter accretes onto horizontal field lines, and showed that inward field advection 
can compete with outward Ohmic diffusion only when the accretion rate is larger than 
10% of the Eddington rate; Payne & Melatos (2004) constructed global hydromagnetic 
equilibrium sequence of a neutron star accreting at the magnetic poles (see also Melatos 
& Phinney 2001); Lovelace et al. (2005) considered a model in which the magnetic field 
is screened by the current in the infalling plasma outside the star. In any case, it is 
not clear that such field burial is effective or permanent, as the field tends to resurface 
due to various instabilities (e.g. buoyancy). If the original field is sustained by currents 
in the crust, then accelerated Ohmic decay due to accretion-induced heating may be a 
viable mechanism (e.g., Urpin & Geppert 1995; Konar & Bhattacharya 1998). 

A different type of accretion-induced field reduction mechanism relies on the 
interaction between superfluid neutron vortices and superconducting magnetic fluxoids 
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in the stellar interior (see Srinivasan et al. 1990). When the star spins up, the fluxoids are 
dragged together with the vortices toward the rotation axis, which leads to a reduction 
of the effective dipole moment even though the local field strength at the magnetic 
poles actually increases (Ruderman 1991a, b). Various observational consequences of 
such vortices migration have been reviewed by Ruderman (2004). 

8. Thermal Radiation from Magnetized Neutron Stars 

8.1. Observational Background 

It has long be recognized that thermal, surface emission from isolated neutron stars 
(NSs) can potentially provide invaluable information on the physical properties and 
evolution of NS (equation of state at super-nuclear densities, superfluidity, cooling 
history, magnetic field, surface composition, different NS populations, etc. Sec, 
e.g., Prakash et al. 2001; Yakovlcv & Pethick 2004 for review). In the last few 
years, considerable observational resources (e.g. Chandra and XMM-Newton) have been 
devoted to such study. For example, the spectra of a number of radio pulsars (e.g., 
PSR B1055-52, B0656+14, Geminga and Vela) have been observed to possess thermal 
components that can be attributed to emission from NS surfaces and/or heated polar 
caps (e.g., Becker & Pavlov 2002). Phase-resolved spectroscopic observations are 
becoming possible, revealing the surface magnetic field geometry and emission radius 
of the pulsar (e.g., Caraveo et al. 2004; De Luca et al. 2005; Jackson & Halpern 2005). 
Chandra has also uncovered a number of compact sources in supernova remnants with 
spectra consistent with thermal emission from NSs (see Pavlov et al. 2003), and useful 
constraints on NS cooling physics have been obtained (e.g., Slane et al. 2002; Yakovlev 
& Pethick 2004). 

Surface X-ray emission has also been detected from a number of soft gamma-ray 
repeaters (SGRs) and anomalous X-ray pulsars (AXPs) — these are thought to be 
magnetars, whose radiation is powered by the decay of superstrong [B ^ 10^^ G) 
magnetic fields (see Thompson & Duncan 1995,1996; Woods & Thompson 2005). Fits 
to the quiescent magnetar spectra with blackbody or with crude atmosphere models 
indicate that the thermal X-rays can be attributed to magnetar surface emission at 
temperatures of (3-7) xlO^ K (see, e.g., Juett et al 2002; Tiengo et al. 2002; Patel et 
al. 2003; Kulkarni et al. 2003; Tiengo et al. 2005). One of the intriguing puzzles is the 
absence of spectral features (such as ion cyclotron line around 1 keV for typical magnetar 
field strengths) in the observed thermal spectra. Clearly, detailed observational and 
theoretical studies of surface emission can potentially reveal much about the physical 
conditions and the nature of magnetars. 

Also of great interest are the seven or so isolated, radio-quiet NSs (so-called "dim 
isolated NSs"; see Haberl 2005). These NSs share the common property that their 
spectra appear to be entirely thermal, indicating that the emission arises directly 
from the NS atmospheres, uncontaminated by magnetospheric processes. Thus they 
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offer tlie best liope for inferring tlie precise values of the temperature, surface gravity, 
gravitational redshift and magnetic field strength. The true nature of these sources, 
however, is unclear at present: they could be young cooling NSs, or NSs kept hot by 
accretion from the ISM, or magnetars and their descendants (e.g., van Kerkwijk & 
Kulkarni 2001, Mori & Ruderman 2003; Haberl 2005; Kaspi et al. 2005). Given their 
interest, these isolated NSs have been intensively studied by deep Chandra and XMM- 
Newton observations. While the brightest of these, RX J1856. 5-3754, has a featureless 
spectrum remarkably well described by a blackbody (Drake et al. 2002; Burwitz et 
al. 2003), absorption lines/features at -E ~ 0.2-2 keV have recently been detected from at 
least four sources, including IE 1207.4-5209 (0.7 and 1.4 keV, possibly also 2.1, 2.8 keV; 
Sanwal et al. 2002; DeLuca et al. 2004; Mori et al. 2005), RX J1308.6+2127 (0.2-0.3 keV; 
Haberl et al. 2003), RX J1605.3+3249 (0.45 keV; van Kerkwijk et al. 2004) and RX 
J0720.4— 3125 (0.27 keV; Haberl et al. 2004), and possibly two additional sources (see 
Zane et al. 2005). The identifications of these features, however, remain uncertain, with 
suggestions ranging from electron/ion cyclotron lines to atomic transitions of H, He or 
mid-Z atoms in a strong magnetic field (Ho & Lai 2004; Pavlov & Bezchastnov 2005; 
Mori et al. 2005). These sources also have different X-ray light curves: e.g., RX J1856.5- 
3754 and RX J1605+3249 show no variabihty (pulse fraction 1-3%); RX J0720-3125 
shows a single-peaked P = 3.39 s pulsation of amplitude ~ 11%, with the spectral 
hardness and line width varying with the pulse phase; RX J 1308-1-2127 has a double- 
peaked pulsation with P = 10.3 s and amplitude ~ 18%. Another puzzle concerns the 
optical emission: For at least four of these sources, the optical counterparts have been 
identified, but the optical flux is larger (by a factor of 4-10) than the extrapolation from 
the black-body fit to the X-ray spectrum (see Haberl 2005). 

The preceding paragraphs highlight the great observational progress made in recent 
years on the study of NS surface emission. These observations can potentially tell us 
much about the physics and astrophysics of NSs. Future X-ray telescopes ( Constellation- 
X and XEUS) will have great capability of studying NS surface emission. In order to 
properly interpret the current and future observations, it is crucial to have a detailed 
understanding of the radiative properties of the outer layers of NSs in the presence of 
intense magnetic fields, and to model the emergent thermal radiation spectra from the 
NS. 

8.2. Modeling Neutron Star Surface Radiation 

Thermal radiation from a magnetized NS is mediated by the thin atmospheric layer (with 
scale height 0.1-10 cm and density ~ 10~^-10^ g cm~^) that covers the stellar surface. 
The physical properties of the atmosphere, such as the chemical composition, equation 
of state, and especially the radiative opacities, directly determine the characteristics of 
the thermal emission. While the surface composition of the NS is unknown, a great 
simplification arises due to the efficient gravitational separation of light and heavy 
elements (see ^6.4p . The strong magnetic field makes the atmospheric plasma anisotropic 
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and birefringent. If the surface temperature is not too high, atoms and molecules may 
form in the atmosphere. Moreover, if the magnetic field is sufficiently strong, the NS 
envelope may transform into a condensed phase with very little gas above it ( §fi.4j) . A 
superstrong magnetic field will also make some quantum electrodynamics (QED) effects 
(e.g., vacuum polarization) important in calculating the surface radiation spectrum (see 
below). 

The first models of zero-field NS atmospheres were constructed by Romani (1987). 
Further works used improved opacity and equation of state data from the OPAL 
project for pure hydrogen, helium and iron compositions (Rajagopal & Romani 1996; 
Zavlin et al. 1996; Gansicke et al. 2002; Pons et al. 2002). These models may be 
applicable to weakly magnetized {B ^ 10® G) NSs. So far most studies of magnetic 
NS atmospheres have focused on hydrogen and moderate field strengths of 5 ~ 10^^- 
10^3 G (e.g., Shibanov et al. 1992; Zane et al. 2000; Zavlin & Pavlov 2002). These 
models take into account the transport of different photon modes through a mostly 
ionized medium. The opacities adopted in the models include free-free transitions and 
electron scattering. However, since a strong magnetic field greatly increases the binding 
energies of bound species (see atoms, molecules, and other bound states may have 
appreciable abundances in the atmosphere (see ^6.4|) . Thus these magnetic atmosphere 
models are expected to be valid only for relatively high temperatures (T ^ a few x 10^ K) 
where hydrogen is almost completely ionized. As the magnetic field increases, we expect 
these models to break down at even higher temperatures as bound atoms, molecules and 
condensate become increasingly important. Models of magnetic iron atmospheres (with 
B ~ 10^^ G) were studied by Rajagopal et al. (1997). Because of the complexity in the 
atomic physics and radiative transport, these Fe models are necessarily crude. 

NS atmospheres with superstrong {B ^ 10^^ G) magnetic field have been studied 
recently (Ho & Lai 2001,2003,2004; Ozel 2001; Zane et al. 2001). Currently most models 
treat ionized H and He atmospheres, and include full angle-dependent transport of the 
photon polarization modes and ion cyclotron resonance in the opacity. The effect of 
vacuum polarization has been studied, and is found to be important for determining 
the thermal spectra and polarizations of magnetized NSs (see ^8.3|) . Recently, NS 
atmosphere models that include self-consistent treatment of the thermodynamics and 
opacities of bound H atoms in strong magnetic fields were constructed (see §8.4j) . 

In order to confront theoretical atmosphere models with observational data 
(spectra, light curves, and in some cases phase-resolved spectra), it is necessary to 
calculate synthetic spectra from the whole stellar surface, taking into account the 
effect of gravitational redshift and light-bending. This can be accomplished using the 
standard, well-developed procedure (Pechenick et al. 1983; Beloborodov 2002). To 
calculate synthetic spectra, one must know the distributions of magnetic field (both 
magnitude and direction) and the effective temperature over the stellar surface; these 
are clearly model-dependent (see Zane & Turolla 2006). For a passively cooling NS 
with a given field geometry, the surface temperature distribution can be obtained from 
the results of NS heat conduction calculations (see ^6.3|) . provided no extra heating 
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Figure 14. A schematic diagram illustrating how vacuum polarization-induced mode 
conversion affects the emergent radiation from a magnetized NS atmosphere. The 
photosphere is defined by where the optical depth (measured from the surface) is 
2/3 and is where the photon decouples from the matter. The left side applies to 
the "normal" field regime [i? < 7 x lO^'^ G; see Eq. H14U|l ]. in which the vacuum 
resonance lies outside the photospheres of the two modes. The right side applies to 
the "superstrong" field regime {B ^ 7 x 10^^ G), in which the vacuum resonance lies 
between the two photospheres. In this diagram, complete adiabatic mode conversion 
is assumed across the vacuum resonance. From Ho & Lai (2004). 



(e.g., associated with particle bombardment in the polar cap) is present. Obviously, any 
spectral feature is expected to be broadened if different parts of the NS surface with 
different field strengths contribute similarly to the observed flux (e.g., Zane et al. 2001; 
Ho & Lai 2004). 

8.3. Effect of Vaccum Polarization on Neutron Star Thermal Spectrum and 
Polarization 

Vacuum polarization can dramatically affect the surface emission from magnetized NSs 
(Lai & Ho 2002,2003a). As discussed in ^4.2M.3[ quantum electrodynamics predicts 
that in strong magnetic fields the vacuum becomes birefringent. In a magnetized 
NS atmosphere, both the plasma and vacuum polarization contribute to the dielectric 
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property of the medium, and a "vacuum resonance" arises when these two contributions 
"compensate" each other. When a photon propagates outward in the NS atmosphere, it 
may convert from one mode into another as it traverses the vacuum resonance. Because 
the two photon modes have very different opacities, the vacuum polarization-induced 
mode conversion can significantly affect radiative transfer in magnetized atmospheres. 
When the vacuum polarization effect is neglected (nv), the decoupling densities of the 
0-mode and X-mode photons (i.e., the densities of their respective photospheres) are 
approximately given by (see Lai & Ho 2002) 

Po,nv ^ 0.42 Tg-^/^Ef '5-1/2 g ^^-3 (138) 
px,nv ^ mT,-'/'El/'S-'/^B,, g cm-3, (139) 

where Tq = T/(10^K) and S* = 1 — e"^/'^^^. Thus the X-mode photons are produced 
in deeper, hotter layers in the atmosphere than the 0-mode photons. When vacuum 
polarization is taken into account, the decoupling densities can be altered depending 
on the location of the vacuum resonance pv relative to po,nv and px.nv For "normal" 
magnetic fields, B < Bi, with 

Bi ^ 6.6 X 10^2 Tq^^^E^^^^S-^/'^ G, (140) 

the vacuum resonance lies outside both photospheres (py < po,nv < Px,nv), and 
the net thermal spectrum is not affected by the vacuum resonance. For superstrong 
magnetic fields, B > Bi, the vacuum resonance lies between these two photospheres 
(Po,nv < Pv < Px,nv), and the effective decoupling depths of the photons are changed. 
Indeed, we see from Fig.Elfhat mode conversion makes the effective decoupling density 
of X-mode photons (which carry the bulk of the thermal energy) smaller, thereby 
depleting the high-energy tail of the spectrum and making the spectrum closer to 
black-body (although the spectrum is still harder than black-body because of nongrey 
opacities)*. This expectation is borne out in the atmosphere model calculations (Ho 
& Lai 2003). Another important effect of vacuum polarization on the spectrum is the 
suppression of proton cyclotron lines and other spectral lines (Ho & Lai 2003; Ho et 
al. 2003). The physical origin for such line suppression is related to the depletion 
of continuum flux, which makes the decouphng depths inside and outside the line 
similar. It was suggested (Ho & Lai 2003) that the absence of cyclotron lines in the 
quiescent spectra of several magnetars (e.g., Juett et al 2002; Patel et al. 2003) is a 
natural consequence of the vacuum polarization effect at work in these systems, and the 
detection of lines in several dim isolated neutron stars (see ^8.1|) is also consistent with 
this picture (Ho & Lai 2004). 

Most studies of radiative transfer in magnetized NS atmospheres rely on solving the 
transfer equations for the specific intensities of the two photon modes (e.g. Meszaros 
1992; Zavlin & Pavlov 2002). These equations cannot properly handle the vacuum- 
induced mode conversion phenomenon. This is because mode conversion intrinsically 

* Note that even when mode conversion is neglected, the X-mode decouphng depth can still be affected 
by vacuum polarization. This is because the X-mode opacity exhibits a spike feature near the resonance, 
and the optical depth across the resonance region can be significant; see Lai & Ho (2002). 
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involves the interference between the modes. In particular, photons with energies 0.3- 
2 keV (this is the energy range in which the bulk of the radiation comes out and spectral 
lines are expected for B ~ 10^^ G) are only partially converted across the vacuum 
resonance. The other problem with the modal description of radiative transport is that it 
is valid only in the limit of large Faraday depolarization (Gnedin & Pavlov 1974), which 
is not always satisfied near the vacuum resonance, especially for superstrong magnetic 
fields (Lai & Ho 2003a). Also, in the presence of dissipation, the two photon modes can 
collapse near the resonance (see §4.3j) . making the modal description meaningless. To 
account for the vacuum resonance effect in a quantitative manner, one must solve the 
transfer equations in terms of the photon intensity matrix (Lai & Ho 2003a) and properly 
take into account of the probability of mode conversion (van Adelsberg & Lai 2006). 
Figure [1^1 depicts an example of such calculation: we see that a superstrong magnetic 
field indeed suppresses the high-energy tail of the thermal spectrum and reduces width 
of spectral lines that may be present. 

We note that even for NSs with "ordinary" field strengths (10^^ < B <7 x 10^^ G), 
vacuum polarization can affect the X-ray polarization signals of the thermal emission 
in a significant way (Lai & Ho 2003b). In particular, the vacuum resonance effect gives 
rise to an unique energy-dependent polarization signature: since the mode conversion 
probability depends on the photon energy f ^4.4|) . the plane of linear polarization at 
E < 1 keV is perpendicular to that at ii^ > 4 keV. By contrast, for superstrong 
field strengths (5 ^ 7 x 10^^ G), the polarization planes at different energies coincide. 
Vacuum polarization is also important for aligning the polarization vectors of photons 
emitted from different patches of the NS, thus ensuring an appreciable net polarization 
fraction at the observer (Heyl et al. 2003; see §4.4|1 . The detection of polarized X-rays 
from neutron stars can provide a direct probe of strong-field quantum electrodynamics 
and constrain the neutron star magnetic field and geometry. 

8.4- Partially Ionized Atmospheres and Radiation from Condensed Surfaces 

For sufficiently strong magnetic field and/or low temperature, bound species have non- 
negligible abundances in the NS atmosphere. Early considerations of partially ionized 
atmospheres (e.g., Rajagopal et al. 1997) relied on oversimplified treatments of atomic 
physics and nonideal plasma effects in strong magnetic fields. In recent years, significant 
progress has been made in studying partially ionized hydrogen plasmas: The binding 
energies and radiative transition rates, incorporating the center-of-mass motion effects 
(see ^3.111 . have been calculated for the H atom (Potekhin 1994; Potekhin & Pavlov 
1997); these atomic data have been implemented in calculations of thermodynamic 
functions (including ionization equilibrium; see ^6.4p and radiative opacities (Potekhin 
et al. 1999; Potekhin & Chabrier 2003,2004). The bound species contribute to the 
bound-bound and bound-free opacities, which, because of the center-of-mass motion 
effect, do not have a narrow width or sharp edge. In addition, the bound species affect 
the dielectric tensor of the medium and hence the polarization properties of the normal 
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Figure 15. The spectrum of a fully ionized hydrogen atmosphere with B = 5x 10^^ G 
and Teff = 5 x 10^ K. The solid line shows the result when the vacuum polarization 
effect is properly taken into account (including partial mode conversion at the vacuum 
resonance), the dotted line shows the result when the vacuum effect is turned off, and 
the dashed line is for a blackbody with T = 5 x 10^ K. The Ed = 0.63 keV ion cyclotron 
feature is prominent in the "no vacuum" curve, but is suppressed in the "partial conv" 
curve. From van Adelsberg & Lai (2006). 

modes; these were studied by Potekhin et al. (2004) using the Kramers-Kronig relation 
between the real and imaginary parts of the plasma polarizability. All these physical 
ingredients have been incorporated into NS atmosphere models to produce the spectrum 
of surface emission (Ho et al. 2003; Potekhin et al. 2004). An example is shown in Fig.ll6| 
clearly demonstrating the importance of the partial-ionization effects. 

As discussed in ^6.4[ for sufficiently low temperatures and/or strong magnetic fields, 
the atmosphere may undergo a phase transition into a condensed state. Thermal 
emission from such a condensed (Fe or H) surface was studied by van Adelsberg 
et al. (2005) (cf. Brinkmann 1980; TuroUa et al. 2004; Perez- Azorin et al. 2005). 
Obviously, the standard radiative transfer equation does not apply inside the condensed 
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Figure 16. Spectral flux as a function of the photon energy E. Left panel: The case 
of a partially ionized hydrogen atmosphere model (solid line) at B — 10^^ G (field 
normal to the surface) and Toff 10^ K is compared with the fully ionized model 
(dashed line) and with the blackbody spectrum (dotted line). Adapted from Ho et 
al. (2003) and Potekhin et al. (2004,2006). Right panel: The cases of condensed Fe 
surface {B = lO^^ G, short-dashed hne; lO^^ G, dotted Hue) and H surface {B = 10^^ G, 
long-dashed line) at temperature T = 10^ K, compared with the blackbody spectrum 
(solid line). Adapted from van Adelsberg et al. (2005). 

phase. Instead, we can apply the Kirchkoff 's law for a macroscopic object (generalized 
for polarized radiation; see §3.1 of van Adelsberg et al. 2005): We calculate the 
(polarization-averaged) reflectivity of the surface for a incident ray with incident angle 6; 
then the emission intensity of the surface is simply 1^ = [1 — Ru{0)] By{T), where By{T) 
is the Planck function. To calculate Ru{d), we can approximate the dielectric tensor 
of the condensed surface using the quasi-electron gas approximation, with and without 
the ion response. For a smooth condensed surface, the overall emission is reduced from 
the blackbody by less than a factor of 2, The spectrum exhibits modest deviation from 
blackbody across a wide energy range, and shows mild absorption features associated 
with the ion cyclotron frequency and the electron plasma frequency in the condensed 
matter. Examples of the spectrum for different models of the surface are shown on the 
right panel of Fig. If the surface is rough (as is likely in the Fe case), the surface 
reflectivity will be close to zero, making emission spectrum even closer to blackbody. It 
is possible that the almost perfect blackbody spectrum in the X-ray band observed from 
RX J1856. 5-3754 (see ^8.1|) is a result of such condensed surface emission, although an 
expanation for the excess optical emission is still called for. 
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8. 5. Resonant Scatterings of Surface Emission 

An interesting issue concerns the possibility that thermal emission from the NS surface 
may be modified by scatterings in the magnetospheric plasma. Since a thermal X-ray 
photon has energy less than the electron cyclotron energy {Ece = 11.6 -B12 keV) at 
the stellar surface, as the photon passes through the magnetosphere, it goes through 
the electron cyclotron resonance {E = Ece) where enhanced scattering occurs. For 
the Goldreich- Julian charge density in the closed magnetosphere, the optical depth of 
such scattering is negligible, thus the surface emission will not be modified (Rajagopal 
& Romani 1997). Ruderman (2003) suggested that a larger density of pair plasmas 
may be present in the magnetosphere, maintained by conversion of 7-rays from the 
NS's polar-cap and/or out-gap accelerators. Partly motivated by the observational 
puzzles associated with the spectra of isolated NSs, he further suggested that resonant 
scattering of thermal photons will result in a Planck-like spectrum that is modified 
from the seed surface spectrum (see also Wang et al. 1998; Lyutikov & Gavriil 
2006). Also, if the plasma is sufficiently hot, it can Compton upscatter the thermal 
photons. While this scenario may be plausible for energetic rotation-powered pulsars, 
the situation is not clear for the dim, radio-quiet isolated NSs, given the absence of 
any nonthermal emission signature from these stars. For active magnetars, a corona 
consisting mainly of relativistic electron-positron pairs can be generated by crustal 
magnetic field twisting/shearing due to starquakes, with the plasma density much higher 
(by a factor of order c/Qr) than the Goldreich- Julian value (Thompson et al. 2002; 
Thompson & Beloborodov 2005). Such a corona plasma can significantly modify the 
surface radiation by multiple scatterings; this may explain the 2 — 10 keV soft tail of 
the magnet ar surface emission. 

9. Non-Thermal Radiation of Rotation-Powered Pulsars 

There is abundant evidence that rotation-powered pulsars must be capable of particle 
acceleration to energies of at least 10 TeV. Pulsed 7-rays above 100 MeV have been 
detected from seven pulsars with high significance and about five others with less 
significance (Kanbach 2002). Some 60 pulsars have been detected at X-ray energies 
with about 30 having detectable pulsations (Kaspi et al. 2005). Both the acceleration 
voltage and the high-energy luminosities required can originate from the rotation of 
a dipole field, if the surface magnetic fields Bg are in fact those required to produce 
the observed rate of spin-down by magnetic dipole torques [see Eqn (jT}]. The vacuum 
potential drop across the open field lines is 



where B12 = G and P is the pulsar spin period (in s). The fraction of this 

potential available for particle acceleration varies among different models and from 
source to source (see Harding 2005 for review), but the overall picture seems energetically 
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self-consistent (although unsolved problems such as global current closure still exist). 
Two accelerator sites have been studied in most detail and two main types of model have 
developed. Polar cap models (Daugherty & Harding 1996, Usov & Melrose 1996) are 
based on particle acceleration beginning near the magnetic poles, where the magnetic 
fields are high and one-photon pair creation ( ^5.4.111 usually dominates over two-photon 
pair creation. Outer gap models (Cheng et al. 1986, Romani 1996) are based on 
acceleration in vacuum gaps in the outer magnetosphere where the magnetic fields are 
too low for one-photon pair creation, but pairs may be created by interactions of 7-rays 
produced by gap-accelerated particles with X-ray photons from either the gap or the 
hot neutron star. 

9.1. Electromagnetic cascades 

9.1.1. Polar cap cascades The combination of high particle energies and magnetic 
field strength approaching, and in some sources exceeding, the critical field allow very 
favorable conditions for pair creation and electromagnetic cascades (Sturrock 1971) that 
can greatly enhance the number of particles in a pulsar magnetosphere (Daugherty & 
Harding 1982). The accelerated particles moving along magnetic field lines with high 
Lorentz factors radiate 7-ray photons at very small angles to the field (^0 ~ I/7)) so 
the one-photon pair production rate for these photons is initially zero. However, as 
they propagate through the curved dipole field, their angle increases until the threshold 
condition, esin6' = 2, is reached and the attenuation coefficient becomes large. Each 
member of the pair will be produced in a Landau state whose maximum principal 
quantum number is nmax = 2e'(e' — 2)/i?' (Daugherty & Harding 1983), where e' = esin^ 
is the photon energy in the frame in which it propagates perpendicular to the local 
magnetic field. From the pair production condition, x = ^B'sin6/2 ~ 0.1, an estimate 
for the maximum pair Landau state is 



When the local field B <0.1Bq, pairs will be created above threshold in highly excited 
Landau states and the excitation level is quite sensitive to field strength, n^ax oc 1/B'^. 
The pairs will decay through emission of synchrotron/cyclotron photons, many of which 
will produce more pairs in excited states. A pair cascade can be sustained in such a way 
through several generations, with pair multiplicities reaching as high as 10^ — 10^ (pairs 
per primary electron) (Hibschmann & Arons 2001, Arendt & Eilek 2002). In high-field 
pulsars, where B ~ O.IBq near the neutron star surface, the pair creation attenuation 
coefficient is high enough for pair creation near threshold, so pairs are produced in very 
low Landau states. In this case, cascade pair multiplicities are lower since the number of 
cyclotron photons drop significantly (Baring & Harding 2001). Bound pair production 
(cf. ^5.5|) also becomes important for B ~ O.IBq, which will further lower the pair 
multiplicity. When B ~ Bq, photon splitting dominates over pair production, as was 
discussed in ^5.61 Further reduction of the cascade pair yield in this regime then depends 
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on whether and at what field strength additional splitting modes are possible. If all three 
modes permitted by QED are operating, then a full photon splitting cascade replaces 
the pair cascade and a complete suppression of pairs occurs (Baring & Harding 2001). 
If only the -L— ^-HH mode permitted in the non-dispersive limit is operating, as is more 
likely, then no suppression of pair yields occurs by photon splitting alone. However, the 
II mode photons will create bound rather than free pairs in a strong field (Shabad & 
Usov 1985, 1986), which at least delays the creation of real pairs until the bound pair is 
dissociated by radiation and/or electric fields (Usov & Melrose 1996) or moves to higher 
altitude where the field is smaller. The pair cascade may then resume, but changes 
to pair yields in this case have not been studied. The screening (see next paragraph) 
of the electric field is also delayed by splitting and creation of bound pairs, effectively 
increasing the accelerating potential. The effect of bound pair creation on polar cap 
acceleration has been studied by Usov & Melrose (1996). 

In the case of young pulsars with high surface fields, the pair cascades are limited 
by screening of the electric field by the pairs as they accelerate in opposite directions 
(Arons & Scharlemann 1979). Such pulsars are capable of accelerating primary particles 
to high enough energies (7 ~ 10'^) that their curvature radiation photons can create 
pairs to initiate the cascade. The curvature-radiation initiated cascades produce high 
multiplicities over a small fraction of a neutron star radius and can screen the E\\ on that 
length scale (Harding & Muslimov 2001). Older pulsars with lower magnetic fields have 
difficulty accelerating particles fast enough for their curvature photons to pair produce 
before the magnetic field drops. In this case, inverse-Compton (ICS) photons that are 
produced by the particles at a lower energy (7 ~ 10"^ — 10^), scattering thermal X-rays 
from the neutron star surface, can create pairs at lower altitude (Zhang & Qiao 1996). 
However, the ICS-initiated cascades have much lower multiplicity, since the emission 
rate decreases as the particles' energy increases (in constrast to the CR emission rate 
which increases as with energy), and are not robust enough to screen the E\\ (Harding 
& Muslimov 2002). These cascades are thus limited by the increasing energy of the 
accelerating particles and the E\\ ultimately saturates at several polar cap radii by 
geometric effects. In the case of millisecond pulsars having extremely low magnetic 
fields, but very short periods, two-photon pair creation through interaction of curvature 
radiation or ICS photons with thermal X-ray photons from the neutron star surface 
can dominate over one-photon pair creation (Zhang & Qiao 1998, Harding et al. 2002). 
Ultimately, as the pulsar ages the polar cap potential drop becomes too low to produce 
pairs of any kind (Harding & Mulsimov 2002, Harding et al. 2002). The resulting 
"death line" in P-P space accounts for the abrupt decrease in numbers of radio pulsars 
at periods longer than a few seconds (see Figure [T]). 

9.1.2. Outer gap cascades Pair cascades in the outer magnetosphere of a pulsar, where 
magnetic field strengths are in the range B ~ 10^ — 10^ G, must be sustained by two- 
photon pair creation. The cascade- initiating particles are accelerated to high Lorentz 
factors in near- vacuum electric fields that develop as a result of the depletion of charge 
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beyond the null-charge surface Q- B — (Cheng, Ho & Ruderman 1986). They radiate 
curvature, or more accurately, synchro-curvature (Zhang & Cheng 1996) 7 rays tangent 
to magnetic field lines within the gap and propagate until they interact with soft photons, 
cither from the hot neutron star surface or from lower energy gap radiation, to create 
a pairs. The pair accelerates in opposite directions, emitting more synchro-curvature 
radiation photons which produce more pairs and the secondary particles supply current 
through the gap. The multiphcity of outer-gap cascades could be quite high (Hirotani, 
priv. comm. 2005), with most of the pairs produced by collisions between downward- 
going gamma-rays and surface X-rays, with magnetic pair creation being important near 
the neutron star surface. These cascades are limited by screening of the E\\ in directions 
both parallel and perpendicular to the magnetic field (Cheng et al. 2000, Hirotani & 
Shibata 2001), taking into account the curvature of the field lines. The gaps can also 
be self-limited if the thermal X rays supplied by the hot neutron star surface require 
heating by gap-accelerated particles flowing downward (Cheng 1994, Zhang & Cheng 
1997). However, neutron star coohng that releases the latent heat of formation, provides 
most of the surface X rays in the youngest pulsars. 

Older pulsars cannot sustain pair cascades in an outer gap, because the E\\ is smaller 
and even the potential drop across the entire outer magnetosphere is not sufficient for 
pair creation. An outer-gap death line in P-P space is thus predicted (Chen & Ruderman 
1993), which may also be a function of pulsar inclination angle (Zhang et al. 2004), and 
suggests that pulsars older than r ~ 10^ yr do not have outer gaps or the associated 
high-energy emission. 



9.2. Observable signatures of high-field physics 



There are several predicted features of non-thermal pulsar emission that would be 
signatures of physical processes operating exclusively in magnetic fields near Bq. The 
spectra of polar cap cascade emission will display high-energy cutoffs due to one-photon 
pair creation attenuation (Daugherty & Harding 1982). Such cutoffs are distinguishable 
by their very sharp, super-exponential shape, so that the observed spectra will have the 
form 

/(e) = Ae-« exp [-Ci, exp {-ell/e)] (143) 

where a is the power-law index, Ci^ = 0.2{a/ X){B'R'^ / p) and e^sc is the escape energy 
from the magnetosphere for photons of energy e (Harding et al. 1997, Harding 2001), 

ell - 2 GeVPV2 g^^^ ^^,|o.i, 5- (^)'| (144) 

or the energy above which photons emitted from a radius ro can escape the 
magnetosphere without pair attenuation. The above estimate of the escape energy 
assumes that the condition for pair attenuation is x ~ 0.2 for B' <0.1 and esin^ = 2, 
the threshold condition, for ~ 0.1, that sin ^ ~ it!/p, where p is the field fine radius of 
curvature at the polar cap rim and that p and B' are evaluated at the photon emission 
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radius r. The cutoff in the spectrum of polar cap cascades that are extended over more 
than a stellar radius will be somewhat more gradual than Eqn ()143|) due to variations 
in both the field strength and radius of curvature. Nevertheless, the cutoff should still 
be distinguishable from a simple exponentional. 

For pulsars with high magnetic fields i?' ~ 1, photon splitting dominates the 
attenuation of polar cap cascade photons in the _L mode (see ^5.6|) . while || mode photons 
are attenuated by pair production. The photon splitting cutoff occurs at roughly the 
escape energy for photon splitting, 



where B' = {ro/R)~^ is the local field strength. At high fields the photon splitting 
escape energy will be lower than that for pair production. The photon splitting cutoff 
will also have a more gradual shape since the attenuation coefficient is power law rather 
than exponential. Because the _L mode photons are attenuated at a lower energy than 
the II mode photons, the spectrum will be 100% polarized in the || mode in the energy 
band between el^^ and el]^ (Harding et al. 1997). This band will grow as the field 
strength increases, as the e*^^ decreases and the field dependence of el]^ saturates for 
5' ~ 1. If such strong polarization in 7-ray pulsar spectra near the high energy cutoff is 
detected by future X-ray and 7-ray polarimeters, it would be a clear signature of photon 
splitting. 

The above estimates for pair creation and photon splitting escape energies assumes 
flat spacetime. Near a neutron star, depending on its compactness (i.e. equation of 
state), general relativistic effects including photon red-shift, light bending and changes 
to the dipole field are important for photon attenuation (Gonthier & Harding 1994). 
For standard neutron star compactness, curved spacetime effects will lower the escape 
energies for both pair creation and photon splitting by a factor of about 2 (Harding et 
al. 1997, Baring & Harding 2001). 

10. Radiation of Magnetars 

SGRs and AXPs were both discovered several decades ago, but they have only recently 
been recognized as similar objects (for detailed review, see Woods & Thompson 2005, 
Heyl 2005). SGRs were first detected around 1979 as 7-ray transients and were thought 
to be a type of classical 7-ray burst. They undergo repeated bursts with several tenths 
of second duration and average energy 10^° — 10'^^ erg, and their bursting often occurs in 
episodes spaced years apart. They more rarely undergo giant superflares of total energy 
_ ]^g47 gj,g^ consisting of an initial spike of duration several hundred ms followed by a 
longer decay of duration several hundred seconds showing pulsations. Such superflares 
have been observed in three SGR sources, SGR0526-66 (the famous 5th March 1979 
event), SGR1900+14 (Hurley et al. 1999) and very recently in SGR1806-20 (Hurley 
et al. 2005; Palmer et al. 2005). In 1998, SGR1806-20 was discovered to have 7.47 s 
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pulsations in its quiescent X-ray emission (Kouveliotou et al. 1998) (and a large P that 
implied a huge surface magnetic field of ~ 10^^ G if due to dipole spin-down. Quiescent 
periodicities of 8 s and 5.16 s and large P were subsequently detected in SGR0526-66 
and SGR1900+14, implying similarly high surface magnetic fields. In all three sources, 
the quiescent periods are the same periods seen in the decay phases of their superflares. 
The quiescent pulse profiles are very broad and undergo dramatic changes just before 
and after superfiares. The profiles are often more complex, with multiple peaks before 
fiares, changing to more simple single peaks profiles following the fiares. All of the SGRs 
lie near the Galactic plane and are thought to have distances around 10-15 kpc (except 
for SGR0526-66, which is in the LMC). 

The quiescent spectra of AXPs and SGRs consist of a thermal component fit by ~ 
0.5-1 keV blackbodies and one or more non-thermal components, with X-ray luminosity 
Iq35 _ ]^g36 gj,g g-i_ "(jj^^j} recently, the non-thermal spectra below 10 keV were fit with 
steep power laws having indices F ~ 2 — 4 (spectral fiux oc e~'"). When INTEGRAL 
and RXTE recently measured the spectra above 10 keV for the first time hard, non- 
thermal components were discovered in three AXPs, and also SGR 1806-20. In two 
of the AXPs, the differential spectra between 10 keV and 50 keV are extremely flat: 
IE 1841-045 (Kuiper et al. 2004) has a power-law index of s = 0.94 and 4U 0142+61 
(Hartog et al. 2004) displays an index of s = 0.45, both much flatter than the steep non- 
thermal components in the classic X-ray band. RXS J1708-40 possesses a slightly steeper 
continuum with s = 1.18. The non-thermal tail of quiescent emission in SGR 1806-20 
is similarly pronounced (Mereghetti et al. 2005, Molkov et al. 2005)), but somewhat 
steeper, with an index of s = 1.6 — 1.9 extending beyond 100 keV. The luminosity in hard 
x-rays (~ 100 keV), ~ 10^^ erg s~^, even exceeds the thermal luminosity from the stellar 
surface. Such hard non-thermal components require continuous particle acceleration 
during the quiescent state. 

Although it was proposed early on (Ramaty et al. 1980, Katz 1982) that SGRs were 
neutron stars with strong magnetic flelds in the range Bs ~ 10^^ — 10^^ G to conflne the 
burst radiation, it was sometime later that surface magnetic flelds exceeding 10 ^'^ G were 
actually predicted for these objects (Duncan & Thompson 1992; Thompson & Duncan 
1995,1996). Duncan & Thompson (DT), exploring the origin of neutron star magnetic 
flelds, conjectured that in dynamo models magnetic flelds 3 x 10^'' G {Pi/1 ms)~^ could be 
generated in nascent neutron stars with initial periods Pi in the ms range (see Thompson 
& Duncan 1993 and gZ21). Paczynski (1992) deduced that SGR0526-66 must have 
a magnetic fleld of 5 x 10^^ G to lower the magnetic scattering cross section f ^5.3|) . 
and therefore raise the Eddington Limit to the luminosity of the repeating bursts at a 
distance of the LMC (see also Miller 1995). This is coincident ally the fleld required if 
its spin-down age at rotation period of 8 s is compatible with the age of its associated 
supernova remnant N49 (Duncan & Thompson 1992). Fields of this magnitude will 
decay of much shorter timescales due to the onset of ambipolar diffusion ([J7j) 
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Figure 17. Evolution of effective neutron star surface temperature when the magnetic 
field decays through the irrotational mode, for different initial magnetic fields. The 
light lines are results for no field decay and the long-short dashed line shows the cooling 
evolution for an unmagnetized neutron star. From Heyl & Kulkarni (1998). 



Diffusion of magnetic ffux out of tlie NS core on tliese timescales provides tlie power to 
magnetars in tlie DT model. Magnetar-strengtli fields also apply higher stresses to the 
stellar crust, so that the yield strain can exceed the crustal strength. This property is 
responsible for the small SGR and AXP bursts in the DT model (Thompson & Duncan 
1996). If a toroidal component of the field -Score > 10^^ G develops in the interior of 
the star, it can twist the external field (Thompson & Duncan 2001). Such action can 
cause the superfiares if the twisted field lines reconnect. The field line twists following 
starquakes (and the associated X-ray burst activities) may last for years and lead to 
significant heating of the near vicinity of the star (Thompson et al. 2002; Thompson 
& Beloborodov 2005). Finally, due to the much faster heat transport in very strong 
magnetic fields, there is a modest increase in the heat flux through the crust Heyl 
& Kulkarni 1998, see Figure ITTjl. All these effects may explain the much hotter surface 
temperatures of magnetars and the high quiescent X-ray emission. 

The magnetar model has had success in accounting for the energetics of the burst 
and quiescent emission in terms of the response of the neutron star to the enormous 
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stresses of field decay. According to the DT model (Thompson & Duncan 1995), the 
magnetar superflares results from reconnection of sheared or twisted external field lines, 
leading to particle acceleration and radiation of hard emission (see also Lyutikov 2006). 
The estimated total energy of such events, 

^«3.4><10-e.g(^)^ (147) 

is similar to observed luminosities of superflares. The smaller bursts result from cracking 
of the crust, which is continually overstressed by diffusion of magnetic flux from the NS 
interior. The shaking of magnetic footpoints then excites Alfven waves that accelerate 
particles. The energy radiated in such events would be 

-i41 / S \ 111 '^max 



^,,„.10"erg^_|_J . (148) 

where / is the length scale of the displacements, Bs is the crustal field and 6'max is the 
yield strain of the crust. The quiescent emission in the DT model is power by magnetic 
field decay through conduction of heat from the core. The NS crust is heated to a 
temperature of 



- crust 



1.3 X lO^iT (^) (149) 
where Tcore is the core temperature, and luminosity 

L.-6xlO^W-^(^)". (150) 
Model atmospheres for the quiescent thermal emission from magnetars was discussed in 

m 

Magnetar models have so far not had complete success in accounting for the 
spectrum of the burst radiation. In the DT model, 7 rays are produced by scattering of 
X-rays from the neutron star surface by particles accelerated after the energy released 
in crust cracking. In the equatorial region, the 7 rays are attenuated by photon splitting 
and pair creation (assuming that only the -L^|||| splitting mode operates), forming a 
pair plasma that is optically thick to Compton scattering at an equilibrium temperature 
of ~ 100 keV (Thompson & Duncan 1995) and confined by the magnetic field. Such an 
estimate results from the deposition of an energy ~ lO'^^ erg into a volume ~ -Rmax; 
with R 

max ~ 10 R. Since the observed burst spectra are fitted with much lower 
temperatures (~ 10 — 20 keV), we cannot be viewing the pair plasma directly, and 
perhaps the final spectrum is formed at small optical depths. Thompson & Duncan 
(1995) argue that if thermal equilibrium of the pair plasma is maintained by photon 
splitting balancing photon merging (its inverse process), then the photon number is 
increased while the temperature decreases to Tgp ~ llkeV, which represents a lower 
limit on the temperature of the emergent SGR burst spectra. But this is a crude estimate 
of the actual temperature of the observed thermal spectrum since scattering and further 
splitting can occur outside the splitting photosphere, the radius at which B ~ Bq. The 
radiative transfer and spectrum of this hot, highly magnetized, and probably dynamic, 
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pair plasma has not been fully modeled. Baring (1995, Harding & Baring 1996) modeled 
spectra resulting from photon splitting cascades assuming that the three splitting modes 
allowed by QED (see §5 .611 operate at magnetar field strengths. The splitting cascade 
spectra are quasi-Maxwellian in shape and the spectral peak saturates in fields above 
~ 5 — IOBq at around 20-30 keV due to the saturation of the photon splitting rate in 
high fields. Although this model could explain the magnetar burst spectral shape, it 
has not been shown that both photon polarization modes can split in magnetar fields. 

Heyl & Hernquist (2005a) propose an alternative model of magnetar radiation that 
appeals to the non-linearity of electrodynamics due to vacuum polarization in a strong 
magnetic field. This model assumes that the power is generated by crustal fractures or 
field re- arrangement, as in the DT model, but postulates that the energy is transported 
primarily by magnetohydrodynamic fast-mode waves rather than by Alfven waves. The 
details of the generation of these waves, whether they are produced directly or through 
interaction of Alfven waves, is not addressed. The fast MHD waves that travel away from 
the neutron star surface with sufficiently large amplitude and short wavelength can break 
down into electron-positron pairs (Heyl & Hernquist 1999) if they form discontinuities 
in the magnetized vacuum on the order of a Compton wavelength. The bulk of the 
pairs form a hot, optically-thick fireball, similar to that produced in the DT model. 
Unfortunately, if the burst spectrum is produced by such an optically-thick pair fireball, 
the actual mechanisms by which it formed are erased and will not be detectable from 
observation. 

Several models have been suggested for explaining the hard component recently 
discovered in quiescent emission. Baring (2004) proposes that resonant Compton 
upscattering of thermal X-rays by accelerated particles in open field region produces 
the quiescent hard emission. Thompson & Beloborodov (2005) have proposed that the 
hard, non-thermal quiescent component is due to the creation of a strong E\\ induced 
by twisting of field in closed region; this leads to synchrotron radiation from electron 
acceleration at high altitude, or free-free emission from a hot (~ 100 keV) transition 
layer between the atmosphere and corona (Thompson & Beloborodov 2005). In the 
model of Heyl & Hernquist (2005b), a small fraction of the energy generated by the 
fast MHD wave breakdown goes into a non-thermal pair component which produces the 
hard quiescent component via a pair-synchrotron cascade. 

10.1. Radio emission and pair suppression in magnetars 

Although transient radio emission has been detected following superflares of 
SGR1900+14 (Frail et al. 1999) and SGR1806-20 (Gaensler et al. 2005), as well as 
the outburst of AXP XTE J1810-197 (Halpern et al. 2005), no persistent, pulsed radio 
emission has been detected from any SGRs or AXPs. The report of radio pulsations 
from SGR1900-I-14 (Shitov et al. 2000) was unconfirmed (Lorimer & Xilouris 2000) and 
probably confused with emission from the nearby radio pulsar PSR J1907-I-0918. It was 
only very recently that radio pulsations were detected from one transient AXP, XTE 
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J1810-197 (Camilo et al. 2006), after its 2003 X-ray outburst. The fact that so few 
magnetars have detected radio pulsations is somewhat of a mystery given the dense pair 
plasmas invoked by magnetar models. It is widely believed that the bright, coherent 
radio emission seen from rotation-powered pulsars is produced by plasma instabilities 
requiring pair plasmas. One must conclude that the pair plasma instabilities operating 
to produce radio emission in rotation-powered pulsars are not effective in magnetar 
fields, that the radio emission is absorbed, or that magnetars do not produce dense pair 
plasmas. The first option is complicated by the fact that radio pulsars with magnetar 
fields exist, so that the inability of pairs to generate the coherent radio emission in 
magnetars must be disrupted by the higher energy densities in their magnetospheres. 
During magnetar bursts, the plasma frequency is almost surely above radio frequencies 
if even a small fraction of the total energy density resides in pairs, so that any radio 
emission would be absorbed. However, in quiescence the particle density is expected 
to be closer to that of the Goldreich- Julian charge of radio pulsars, which gives plasma 
frequencies of 1 - 10 GHz near the surface and falls of r~'^/^ above the surface. Even 
if the particle densities are several hundred times the Goldreich- Julian density, radio 
emission would only be absorbed near the stellar surface. Another reason which has 
been given for the absence of radio pulses is that their small polar caps predict very 
small radio beams that are hard to detect (Thompson & Duncan 1996). But as the 
number of known magnetars increases, we would expect to be in view of a few of the 
radio beams, and a number of radio pulsars with similarly small predicted radio beams 
have been discovered. 

The suppression of pair creation is also a possible explanation for the absence of 
radio emission in magnetars. It was noted that the field strength B ~ Bq at which 
photon splitting begins to dominate the attenuation of photons near pulsar polar caps 
occurs at around the dividing line between radio pulsars and magnetars (Baring & 
Harding 1998) and would impose a death line for pair creation at high fields. However, 
the complete suppression of pair creation requires that both photon polarization modes 
undergo splitting (Baring & Harding 2001), but only the _L mode splitting is allowed 
under Adler's (1971) kinematic selection rules f ^5.6|) . Usov (2002) suggested that pair 
suppression could still occur if bound pair creation were also taken into account. Photons 
in the || mode will create bound pairs with both particles in the ground state, which 
are stable against annihilation f ^5.5|) . Photons in the _L mode would create bound pairs 
with one or more particles in an excited state, which can be dissociated by ionization in 
spin-flip transitions, but these photons would split before creating the unstable bound 
pairs. 

11. Accreting X-ray Pulsars 

X-ray pulsars are believed to be strongly-magnetized accreting neutron stars powered 
by kinetic energy of infalling matter from a companion star. The magnetic fields of 
10^^ — 10^^ G are inferred from the pulsations that require anisotropic infall and radiation. 
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the cyclotron lines observed in a number of pulsar spectra (Nagase 1989, dal Fiume et 
al. 2000), and the observed changes in pulse period (Ghosh & Lamb 1979, Bildsten et 
al. 1997). We can get a fairly good estimate of the effective temperature observed in 
X-ray pulsars from T^g — {L^/aAca.pY^'^ — 10 keV, where is the X-ray luminosity and 
Acap is the heated polar cap area. But clearly, since pulsar spectra are not blackbody, 
we must have an accurate description of the radiating plasma in order to model the 
emission. 

The conditions in the emission region are closely tied to the dynamics of the 
accretion flow and infall, and there are two major regimes that depend upon whether or 
not the radiation pressure of the emitting plasma is capable of decelerating the accretion 
flow. The critical luminosity L^rit where this division occurs is not accurately determined 
but can be estimated by requiring that the outgoing radiation pressure balance the ram 
pressure of the infalling matter: 

^crit = — 5^ ~ 10^" ergs \ (151) 

(J|| ii^ 

when a\\ ~ ax and Acap = O.li?^, where a\\ is the magnetic Thomson scattering cross 
section parallel to the field, averaged over the spectrum, polarization and angles of the 
radiation. The value of Lent is about one hundredth of the Eddington limit for spherical 
accretion, and is model dependent since the averaged cross section should include the 
radiation produced throughout the atmosphere in a self-consistent model. The high 
luminosity X-ray pulsars where > L^rit may radiate several orders of magnitude 
above Lcrit and even above the spherical Eddington limit. In low- luminosity sources 
where -C i^crit, radiation pressure is not important in decelerating the accretion flow, 
in which case it may be decelerated by a coUisionless shock above the neutron star 
surface or by Coulomb and nuclear collisions with atmospheric plasma near the surface. 
Most (about two thirds) of the observed X-ray pulsars fall into the high-luminosity 
category. 

The most important effect of the neutron star magnetic field on the processes that 
decelerate the accretion flow and that produce the observed radiation is the quantization 
of particle momentum perpendicular to the field (see 0. In sources where -C Lent, 
accretion flow can convert its kinetic energy to heat either at a coUisionless shock 
front (Basko & Sunyaev 1976) or by Coulomb collisions with thermal electrons and 
protons in the atmosphere (Basko & Sunyaev 1975, Pavlov & Yakovlev 1976). Self- 
consistent models of Coulomb-decelerated acceleration (Harding et al. 1984, Miller et 
al. 1989) result in thin slab-like atmospheres, since Coulomb stopping lengths are only 
~ 50gcm"^ (Kirk & Galloway 1981, Pakey 1990). Models for sources with L^ > Lent, 
where radiation pressure dominates the deceleration, must solve the coupled radiation- 
hydrodynamic flow to find the structure of the atmosphere. In such models, a radiative 
shock stands off several stellar radii above the neutron star surface and the atmosphere 
has more of a cylindrical geometry, as radiation escapes from the sides of a decelerating 
mound of gas below the shock (Langer & Rappaport 1982, Kirk 1985). 

The dominant cooling process in X-ray pulsar atmospheres is resonant 
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Bremsstrahlung (or cyclotron cooling), in which an electron in the ground state is 
collisionally excited and deexcites through spontaneous emission, producing a cyclotron 
photon. The inverse process of resonant free-free absorption, where an electron is excited 
by absorption of a cyclotron photon and collisionally deexcited, can in some cases 
provide significant heating to the atmosphere. So while the vast majority of electron 
excitation-deexcitation events are resonant Compton scattering at the cyclotron energy 
the collisions, although less frequent, control the net production and destruction of 
cyclotron photons in these atmospheres. Another important heating and cooling process 
is Compton scattering, which is resonant at the cyclotron harmonics in a magnetic field. 
While scattering does not produce nearly as much heating and cooling as Bremsstrahlung 
since it does not create and destroy cyclotron photons, it exchanges energy between 
photons and electrons in the atmosphere (Meszaros & Ventura 1979). The presence of 
a strong magnetic field introduces substantial changes in the stopping length due to 
Coulomb scattering. A accreting proton moving along the magnetic field experiences 
a much smaller drag force than in the non-magnetic case because the electrons are 
limited to momentum transfers along the field unless they can be excited to higher 
states. The infalling protons must diffuse in momentum space through many small- 
angle collisions, gradually veering away from the magnetic field direction where their 
drag force will increase and they can decelerate faster. As a result, the stopping lengths 
are very dependent on the proton momentum diffusion and significantly larger than in 
an un- magnetized plasma (see, e.g.. Nelson, Salpeter & Wasserman 1993). 

11.1. Cyclotron lines 

The cyclotron lines seen in the spectra of accreting X-ray pulsars are formed by 
resonant scattering ( %.'S\i of photons with electrons that occupy discrete Landau 
states. If the plasma temperature is less than Tb = es/k, then the electrons occupy 
primarily the ground Landau state [regime (i) of ^ and have a one- dimensional 
thermal distribution of momenta along the magnetic field. At magnetic field strengths 
typical of neutron stars, the cyclotron decay rates from excited Landau states (n > 0), 
r„ = 4 X lO^^ni^fg s~^, are much greater than the rate of coUisional excitation. Unless 
there are other processes that can populate the excited states on faster timescales, 
the bulk of the electrons are expected to occupy the ground state in a one- dimensional 
distribution. Under these conditions, absorption of a cyclotron photon is always followed 
by the emission of another photon, so resonant scattering is more important than 
absorption. As discussed in §5.31 , the Compton scattering cross section has sharp 
maxima at the cyclotron resonance (Eqn [7^) and its harmonics, so that photons 
scatter preferentially at these energies. Since the scattering cross section has a strong 
dependence on angle to the magnetic field, the cyclotron resonance scattering features 
(CRSFs) provide excellent diagnostics of the emission geometry as well the physical 
conditions (temperature and density) of the radiating plasma. 

The formation of CRSFs at the fundamental and at the harmonics proceeds by 
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very different routes. When a photon at the fundamental resonance frequency (Eqn 
[73] with = 1) scatters with an electron in the ground state the scattered photon 
also appears at the fundamental. The scattered photon may have a different angle, 
but unless the continuum photons are beamed, resonant photons from other angles 
scatter into view. Absorption-like features can only form in the fundamental by repeated 
scattering, causing diffusion of photons in energy and angle through the combined effects 
of electron recoil and the natural line width of the resonance (Wasserman & Salpeter 
1980, Wang et al. 1989). Resonant photons can then diffuse from the line core into the 
line wings where the cross section drops to near cr-r and they can then escape, primarily 
in the red wing. When a photon at the first harmonic (A^ = 2) scatters with an electron 
from the same distribution, the scattered photon will appear most of the time at the 
fundamental in fields B' < 0.2, exciting the electron to the n = 1 state. The electron then 
decays through cyclotron emission, spawning an additional photon at the fundamental. 
Although this is technically still a resonant scattering (n = — 2 — 1 followed by 
n = 1 —>■ 0), analogous to Raman scattering in atomic physics, incident photons are 
destroyed in the first harmonic and the line formation here is well approximated by 
absorption (Fenimore et al. 1988, Wang et al. 1988, Harding & Daugherty 1991). 
Thus, the second and higher harmonics are formed in a similar way, with photons 
in the A^th resonance exciting the electron from the ground state to the n — 1 state, 
with the subsequent spawning of n — 1 additional photons at the fundamental through 
single harmonic number transitions. Photon spawning partially fills in the fundamental 
(Wang et al. 1988), so that the line features at the cyclotron fundamental and higher 
harmonics will not appear in the same ratios in the observed spectrum as they do in 
the cross section. In magnetic fields B' ~ 0.2, the dominance of cyclotron decays to the 
ground state f ^5.2|) reduce the rate of photon spawning in the fundamental (Araya & 
Harding 1999) since scattering in the A^th harmonic is more likely to produce a scattered 
photon in the A^th resonance {n = ^ N ^ 0). Formation of the harmonics in these 
high fields therefore proceeds more by diffusion out of the line core, as the fundamental 
is formed, rather than by excitation of the electron. 

As in most line-formation scenarios, the broadening of the CRSFs results from a 
combination of natural line width, which gives a Lorentzian profile, and the Doppler 
effect, which gives a Gaussian profile. Due the one-dimensional momentum distribution 
of the electrons, the line broadening is asymmetric and dependent on viewing angle to 
the magnetic field, 6. First order (non-relativistic) Doppler broadening alone gives a 
line width that is approximately 

Aen ^ eN{2TY^^cose, (152) 

where T is the plasma temperature in units of mc^ and is the resonant energy (Eqn 
[73]). At small viewing angles, the lines are broadened by the full Doppler motion, while 
in the limit of large viewing angles, only transverse (second-order) Doppler broadening 
contributes to the line width. Since the transverse Doppler broadening is a red shift 
(from time dilation), the line profiles become increasingly asymmetric at large incidence 
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Figure 18. Relativistic cyclotron scattering line profiles for B' = 0.24 and thermal 
electron temperature of 31 keV, for different viewing angles, fj, = cos 9. 



angles. The natural line width, due to the finite lifetime of excited Landau states, 
is approximately equal to the cyclotron decay rate from state N which, in the non- 
relativistic limit, is 

Asn ^ ^aB'^N (153) 

At large angles, Aejv provides the main contribution to the broadening on the blue 
side of the line. In a relativistic treatment of line broadening, the cyclotron decay rate 
is dependent on the electron spin, with the rates of spin-flip transitions being smaller 
(Herold et al. 1982). 

Techniques for modeling CRSFs in XRP atmospheres have included either solving 
the equations of radiative transfer using a Feautrier method (Meszaros & Nagel 1985, 
Nagel 1981, Alexander & Meszaros 1991b, Burnard et al. 1988) or Monte Carlo 
simulation (Wang et al. 1988, Araya & Harding 1999, 2000). Given the complexities 
of the fine transfer, these models usually assume a static atmosphere and simplified 
geometries, in particular homogeneous slabs, with the magnetic field perpendicular to 
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the slab normal, or cylinders parallel to the field. An exception is the work of Arons et al. 
(1987) and Burnard et al. (1988) who treat a time-dependent accretion column. Such 
studies have generally confirmed that the shape and relative depths of the fundamental 
and higher harmonics of CRSF are very dependent on viewing geometry as well as 
geometry of the scattering plasma. In particular, CRSFs seen at small angles to the 
magnetic field tend to be broader with weak harmonics, whereas CRSFs seen at large 
angles are narrower, more asymmetric and have stronger features at higher harmonics. 
In slab geometry the depth of the fundamental is largest at large angles to the field, 
where the optical depth is largest, but in cylinder geometry the fundamental feature 
is largest at small angles (Araya & Harding 1999). In many cases the fundamental is 
fianked by emission wings, from photons escaping the line core, and may even appear 
completely in emission, at large angles in cylinder geometry, due to photons diffusing 
from smaller angles. 

Given the sensitive dependence of CRSF characteristics with angle to the field, 
one would expect a great deal of variation in the features with rotational phase of the 
neutron star. This is in fact what is observed in phase resolved spectra of accreting 
XRPs (Heindl et al. 2004, Orlandini & dal Fiume 2001). The line depths as well as 
their energies change with phase, but the harmonic relationships are preserved, as would 
be expected of CRSFs. 

12. Conclusion 

In this article we have reviewed the physics that applies in extremely strong 
magnetic fields, the properties of strongly magnetized neutron stars, atmospheres 
and magnetospheres, and the present status of our understanding and models of the 
astrophysical sources that are their manifestation. While the use of non-relativistic 
quantum mechanics is adequate (even for B > Bq) for describing bound states or any 
processes where the electron stays in the ground Landau level, one must use relativistic 
quantum mechanics (QED) in computing the rates and cross sections for many free 
particle and photon processes. Nearly all of the cross sections for the first and second- 
order free particle processes in strong magnetic fields, such as cyclotron radiation and 
absorption, one photon and two-photon pair creation and annihilation and Compton 
scattering, have been calculated and studied. Only a few of the higher-order processes, 
such as Bremsstrahlung and photon splitting, have been investigated since they become 
increasingly complex. The unusual and interesting properties of photon propagation 
in strong magnetic fields has been a topic of intensive study in recent years, especially 
after the discovery of magnetars. Among all the processes that have been investigated, 
vacuum polarization has been found to be of particular importance. 

A variety of astrophysical sources, that include rotation-powered and accretion- 
powered neutron stars and the magnetars, SGRs and AXPs, are believed to be very 
strongly magnetized neutron stars. The processes that operate only in such strong fields 
may be fundamental to the functioning of these sources. For rotation-powered pulsars. 
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one-photon pair creation is thought to be the primary process attenuating high-energy 
photons, and the created pairs may be critical in the production of the observed radio 
pulsations. We do not yet fully understand though how and why the radio emission 
seems to turn off well before the pulsar spin-down has ceased. Theoretical models show 
that pair production decreases as the pulsar ages, but does it turn off suddenly, or 
are there some threshold properties (e.g. multiplicity or energy spectrum) for radio 
emission? In the case of the magnetars, power levels far exceeding their spin-down 
power requires fast magnetic field decay on timescales only possible through processes 
such as ambipolar diffusion, that become effective in fields above 10^^ G. We do not yet 
understand, however, why these sources are radio quiet. Is pair production suppressed 
by the strong field or are the collective plasma processes disrupted? Theoretical studies 
have so far not been able to find a convincing mechanism for the suppression of the radio 
emission in magnetars. Even more puzzling is the growing number of magnetars and 
radio pulsars having very similar spin properties, and therefore similar implied surface 
magnetic field strengths. Aside from their radio emission properties, radio pulsars 
and magnetars also have very different X-ray emission levels, transient emission and 
glitching behavior. What are the hidden characteristics of these sources that distinguish 
them? Also of interest are the emerging population of "dim" isolated neutron stars 
with apparently pure thermal emission. The spectral features detected in some of these 
sources are exciting, but their identifications remain unclear. The nature of these sources 
is unknown. Could they be descendants of magnetars? 

To find answers to these questions, new ideas and more theoretical work are surely 
needed. But new and more sensitive instruments are on the horizon that will also provide 
some clues as well. The ALFA pulsar survey (Cordes et al. 2005) began operation last 
year at Arecibo, one of the world's most powerful radio telescopes, and this survey is 
expected to discover at least 1000 new radio pulsars, nearly doubling the current number. 
Among the newly discovered pulsars will be many more radio pulsars with magnetar field 
strengths. The Gamma- Ray Large Area Telescope (GLAST) (McEnery et al. 2004), 
due to launch in 2007, will have a point source detection threshold about 20 times below 
that of EGRET and will discover hundreds of new 7-ray pulsars. GLAST will also have 
sensitivity up to 200 GeV and will be able to make very sensitive measurements of 
spectral high-energy cutoffs. Third-generation ground-based Air Cherenkov detectors, 
such as H.E.S.S. (Hinton et al. 2004) in Namibia, have begun operation. They are 
sensitive to 7-rays in the range 50 GeV to 50 TeV and may discover or put important 
constraints on pulsar and magnetar spectra and their nebulae. Further into the future 
arc planned X-ray telescopes, such as Constellation X, XEUS, and X-ray polarimeters, 
such as AXP, POGO and ACT. Polarimeters in particular will be extremely important in 
looking for some of the signatures of very strong magnetic fields that have been discussed 
in this article, such as the vacuum polarization resonance and photon splitting cutoffs. 
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